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Section  1 


INTRODUCTION 


1.1  BACKGROUND 

The  load-bearing  capacity  of  most  alloys  can  be  degraded  by  a 
corrosive  environment.  Two  common  forms  of  environment-aenaitive 
mechanical  property  degradation  are  stress  corrosion  cracking  (SCC) 
and  corrosion  fatigue  (CF).  Both  forms  greatly  influence  the  safety, 
reliability,  and  economics  of  many  components  used  in  military  and 
commercial  power  generation,  transportation,  and  energy  conversion 
systems.  Therefore,  substantial  effort  ha a  been  made  to  understand 
environment-assisted  crack  growth,  to  identify  and  develop  siaterials 
that  have  better  environment  resistance,  and  to  formulate  predictive 
models.1*1^ 

These  efforts  have  resulted  in  various  quantitative  models  being 
proposed  for  predicting  SCC,  CF,  or  other  types  of  environment-assisted 
crack  growth.  Wei  and  Landes, for  examp'*,  considered  the  crack  growth 
rate  in  a  CF  situation  to  consist  of  the  sum  of  the  crack  growth  rate  in 
an  inert  environment  and  the  rate  in  the  environment  of  interest.  The 
Wei  and  Landea  theory  has  been  shown  to  explain  adequately  many  features 
of  the  CF  behavior  of  Ti  and  T1-6A1-4V  in  NaCl  solutions. ^  Furthermore, 
this  work  emphasized  the  need  for  a  better  understanding  of  the  processes 
that  occur  at  the  crack  tip,  particularly  at  threshold  load  levels  for 
crack  growth  and  at  low  losding  frequencies  (<  1  Hz). 

For  SCC  alone,  Logan7  proposed  a  model  for  cracking  in  aqueous 
environments  based  on  a  slip-dissolution  mechanism  at  the  crack  tip. 

This  mechanism,  now  widely  accepted,^-7  relates  crack  growth  to  the 
enhanced  anodic  dissolution  current  at  the  crack  tip  resulting  from  the 
Intermittent  mechanicsl  rupture  of  an  environmentally  created  film  at  the 


Logan *0  model  has  bean  modified  Co  Include  a  variety  of  rate¬ 
controlling  processes  that  Include  repasslvaton  rate,  mlcrocraep, 
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solution  renewal,  and  critical  charge  concepts*  If  the  processes 
occurring  within  the  crack  were  better  understood,  It  would  be  possible 
to  datermlne  the  appropriate  slip-dissolution  model  for  any  given 
material,  environment,  and  load  application* 

Quantitative  modeling  of  the  cracking  process  is  difficult  because 
It  Involves  complex  lntersctions  between  the  cracking  material,  the 
environment,  and  the  applied  load.  Because  specimen  microstructure, 
surface  chemistry,  electrochemistry,  and  loading  conditions  are  expected 
to  Influence  cracking  behavior,  they  should  be  accounted  for  In  attempts 
to  model  crack  growth.  Clearly,  slmultaneoua  characterisation  of  the 
mechanical  and  electrochemical  state  at  the  crack  tip  Is  required  for 
development  of  a  reliable,  predictive  crack  growth  model. 

The  mechanical  state  at  the  crack  tip  la  usually  described  In  these 
models  using  classical  fracture  mechanics  parameters  such  as  the  linear 
elastic  stress  Intensity  (K  parameter)1^* 17  or,  more  recently,  tha 

IB  1  Q 

elastic-plastic  stress  Intensity  (J  para--*ter).  * Unfortunately,  the 
electrochemical  state  of  the  crack  tip  Is  not  so  easily  character¬ 
ized  because,  until  recently,  no  reliable  experimental  method  existed 
that  could  either  provide  measurements  directly  within  the  crack  or 
convert  measurements  made  external  to  the  crack  Into  crack  tip 
parameters.  Some  attempts  have  been  made  to  measure  the  local 
chemistries  that  develop  at  crack  tips  by  using  mlcroanalytlcsl  chemical 
techniques  on  specimens  that  are  actually  cracking  or  on  artlflcal 
crevices  that  had  been  constructed  to  simulate  stress  corrosion  cracks. 
Howevar,  most  of  the  past  crack  growth  models  have  used  electrochemical 
data  obtained  on  bulk  specimens  and  have  assumed  that  such  data,  for 
example,  the  free  corrosion  potential,  are  representative  of  values 
inside  the  crack,  or  these  bulk  values  were  extrapolated  to  crack  tip 
values  by  means  of  some  postulated  model.  However,  crack  growth  rates 

predicted  by  such  models  can  be  several  orders  of  magnitude  different 
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>  fron  those  actually  observed. 
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One  possible  method  for  evaluating  the  electrochemical  state  of 
material  near  a  crack  tip  Is  the  AC  Impedance  technique*  This  technique, 
under  development  for  the  past  several  years, 22-26  ha8  oniy  recently  been 
applied  to  the  problem  of  environmentally  assisted  cracking*  This 
approach  shows  promise  for  determining  crack  tip  electrochemical 
parameters  by  applying  the  transmission  line  model  to  deconvolve  the 
measured  AC  Impedance  spectra  of  growing  cracks* 

1.2  PROGRAM  OBJECTIVES 

The  above  background  discussion  highlights  the  need  to  characterize 
simultaneously  the  electrochemical  and  mechanical  states  In  s  growing 
stress  corrosion  or  corrosion  fatigue  crack*  Thus,  the  objective  of  this 
research  program  was  to  generate  electrochemical  and  mechanical  data  for 
actively  growing  corrosion  fatigue  cracks  and  then  use  these  data  to 
develop  a  preliminary  predictive  model  for  SCC  and  CF* 

To  accomplish  these  objectives  we  used  standard  fracture  mechanics 
methods  and  analyses  to  characterize  the  mechanical  state  of  the  crack  In 
T1-6A1-4V  and  HY-80  steel  In  3.5X  NaCl  solution  and  the  AC  Impedance 
technique  to  characterize  the  electrochemical  state*  The  electrochemical 
and  mechmlcal  data  were  used  to  calibrate  a  transmission  line  model,  and 
this  model  was  then  used  to  predict  quantitatively  the  electrochemical 
and  mechanical  state-at  the  tip  of  an  actively  growing  crack. 

This  report  describes  the  results  of  our  efforts  to  meet  these 
objectives*  Section  2  describes  the  development  of  the  transmission  line 
model*  Section  3  reviews  the  electrochemical  and  mechanical  experimental 
techniques  employed  In  this  program  and  the  results  obtained*  In  Section 
4  the  transmission  line  model  Is  presented  along  with  the  experimental 
data  used  to  develop  a  predictive  model  for  crack  growth.  Section  5 
outlines  the  conclusions  and  recommendations  for  additional  work. 
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Section  2 


DEVELOPMENT  OF  THE  DISCRETE  TRANSMISSION  LINE  MODEL 

In  the  discrete  transmission  line  model  the  stress  corrosion  or 
corrosion-fatigue  crack  Is  treated  as  a  cracked  electrode.  The  solid  and 
electrolyte  phases  In  the  cracked  electrode  are  represented  by  resistance 
elements  parallel  to  the  direction  of  the  corrosion  current  flux,  and  the 
lnterfaclal  resistance  Is  represented  by  an  Impedance  distributed  between 
the  other  two  resistances.  The  electrode  may  be  treated  as  "one- 
dimenslonal“  If  the  potential  (E),  current  flux  (I),  and  reactant 
concentrations  (c)  vary  only  with  the  depth  within  the  crack  and  not  with 
the  lateral  position  along  the  crack.  In  thla  case,  the  local  values  of 
E,  I,  and  c  may  be  replaced  by  their  average  values  In  a  line 
perpendicular  to  the  crack  plane. 

n  7 

Ksenzheck  and  Stender  discussed  the  conditions  that  are  applicable 
to  transforming  a  three-dlmenalonal  problem  Into  a  one-dlmenslonal  one. 
The  mathematical  model  for  the  electrical  representation  shown  In  Figure 
1  requires  solutions  to  a  set  of  differential  equatlona  of  the  following 
form : 
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where  subscripts  e  and  s  refer  to  the  electrolyte  and  solid  phases,  Z  Is 

the  lnterfaclal  Impedance,  x  Is  distance  Into  the  crack,  r  la  the 

resistance,  I  ■  +  I. ,  and  AE  ■  E  -  E  • 

*  8  e»  a  e 

For  the  simplest  set  of  boundary  conditiona  shown  In  Figure  1, 

corresponding  to  an  electrolyte-filled  right  cylindrical  crack  In  which 

rs  and  re  are  Independent  of  x  and  the  solid  and  electrolyte  phase 

current  contacts  are  diametrically  opposed,  the  transmission  line 
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Impedance  Ztl  Is  given  by 


4r„  r  +  (r„2  +  *  2)  (eYl  +  e-*1)  r  r  1 


Jtl 


(r 


re)  (eYl  -  e^1) 


r  +  r 
s  e 


(A) 


where  1  Is  the  current  crack  length,  and 

Y  -  [(r8  +  re)/Zt]1/2  (5) 

The  lnterfaclal  Impedance,  Z^,  In  equation  (5)  Is  equivalent  to  that  at  a 
plane-parallel  electrode.  That  is,  It  contains  the  Information  that 
would  be  available  If  It  were  possible  to  determine  the  electrochemical 
kinetic  properties  of  each  Infinitesimal  element,  dx,  of  the  electrode 
thickness  In  the  x  direction. 
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The  conventional*4  *  equivalent  circuit  for  Z^  Is  shown  In  Figure  2. 
The  components  of  this  circuit  contain  the  electrochemical  kinetic 
parameters  that  are  important  In  evaluating  the  behavior  of  a  cracked 
electrode.  Brief  descriptions  of  these  parameters  follow. 

Double-Layer  Capacitance,  C.  The  double-layer  capacitance,  C,  Is 
directly  proportional  to  the  area  of  the  electrolyte/solld  phase 
Interface.  Thus,  given  a  knowledge  of  the  double-layer  capacitance  per 
unit  area,  C',  under  the  prevailing  conditions  of  reactant  concentra¬ 
tion,  temperature,  and  potential,  one  can  determine  C  by  deconvolving 

the  measured  transmission  line  Impedance.  This  value  of  C  can  then  be 
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used  to  estimate  the  total  wetted  ares  per  unit  electrode  thickness, 
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Reaction  Warburg 


Capacitance,  C 

JA-350 583-29 

FIGURE  2  INTERFACIAL  IMPEDANCE  ELEMENT; 
RANDLES'  EQUIVALENT  CIRCUIT 


which  we  define  as  A'.  Values  of  C'  may  be  available  from  the 
literature  or  can  be  measured  directly  for  a  plane-parallel  electrode  of 
known  area. 


Warburg  Coefficient,  a.  At  moderate  frequencies  (>  1  1L,,  the 

00  0 /  00 

Warburg  Impedance,  can  be  expressed  aa  ’ 


Zw  -  o(l  -  j)w 


1/2 


(6) 


where 


RT 

2  i  - T7T 

n  F  cA"(2D )L/£ 


(7) 


in  which  c  is  the  reactant  concentration  and  A"  expresses  the  electrode 
electrolyte  area  per  unit  length  that  is  actively  engaged  in  reaction. 
The  other  parameters  in  equation  (7)  have  their  usual  meanings.  The 
relationship  between  A'  and  A"  can  be  used  to  define  an  area  efficiency 
that  accounta  for  the  presence  of  passive  zones,  which  do  not  contribute 
to  the  cell  current  during  crack  growth. 

Diffusion  layer  Thickness,  6.  At  the  limit  of  low  frequency,  an 

_  —  29-12 

apparent  reaistive  ahunt  to  the  Warburg  impedance  is  observed.  * 

This  resistance  defines  the  steady-state  DC  operstlng  behavior,  and  its 
value  may  be  used  in  conjunction  with  a  to  calculate  the  effective 
diffusion  layer  thickness,  6.  In  the  derivation  of  equation  (6),  semi¬ 
infinite  linear  diffusion  was  assumed.  However,  in  the  wide  regions 
near  the  crack  mouth,  natural  convection  will  limit  the  thickness  of  the 
diffuaion  layer,  whereaa  in  the  narrower  regions  near  the  crack  tip,  the 

dimension  of  the  diffusion  layer  ia  proportional  to  the  crack  opening. 

jj  10— IS 

Thus,  we  can  account  for  the  steady-atate  DC  behavior  »  J  in  the 
model  of  the  cracked  electrode  by  introducing  the  Nernst  diffusion  layer 
thicknesa  as  a  boundary  condition. 


Interfacisl  Reaction  Resistance,  Rr»  The  reaction  resistance 
together  with  A"  and  6  can  be  used  to  obtain  information  regarding  the 


electrode  kinetics1^*2^*^  within  the  cracked  electrode.  The  exchange 
current  density,  1Q,  can  be  determined  from  the  relatlonahip 


\  -  RT/nFA ' 1Q 


(8) 
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When  U8ed  In  conjunction  with  the  Stem-Geary  relatlonahip,  Rr 
provides  a  measure  of  the  Instantaneous  corrosion  rate  within  the  crack. 

Electrolyte  Resistance,  rg.  The  electrolyte  resistance  per  unit 
length  In  the  cracked  electrode,  rg,  gives  a  measure  of  the  amount, 
composition,  and  distribution  of  electrolyte  In  the  crack.  The 
electrolyte  resistance  Is  controlled  primarily  by  the  crack  width  and  to 
a  smaller  extent  by  the  diffusion  of  reactant  species  from  the  bulk 
electrolyte  and  the  diffusion  of  reaction  products  from  the  crack  tip. 


Solid  Electrode  Resistance,  rn.  The  solid  phase  resistance  per 
unit  length  In  the  cracked  electrode,  rg,  provides  a  measure  of  the 
resistive  component  In  the  cracked  electrode.  Reduction  of  the 
Impedance  data  to  obtain  the  parameters  of  Interest  for  cracked 
electrodes  Involves  applying  nonlinear  regression  analysis  to  sn 
impedance  model  or  equivalent  circuit.  Broers  and  coworkers^®”^® 
demonstrated  the  applicability  of  the  transmission  line  model  for 
characterizing  porous  electrodes  In  molten  carbonate  fuel  cells,  and  we 
believe  that  these  methods  can  be  applied  to  a  growing  SC  or  CF  crack. 

An  additional  impedance  element,  Zt^  ,  has  been  included  In  Figure  1 

to  account  for  the  special  electrochemical  characteristics  of  the  crack 

tip.  For  example,  If  the  crack  grows  by  a  mechanism  of  film  rupture  and 

anodic  dissolution  at  the  crack  tip,  then  Zt^  can  be  represented  as  s 

special  case  of  the  equivalent  circuit  in  Figure  2,  with  R^.  determined  by 

the  exchange  current  density  for  metal  dlsaolutlon  (possibly  involving 

terms  for  precipitation  of  oxide  and  repasalvatlon) .  Under  the 

conditions  of  a  cyclic  mechanical  load,  the  area  term  A'  In  equation  (8) 

« 

must  also  be  given  as  a  periodic  function  about  a  mean  value,  A  .  In 

o 

the  slmpleat  case.  A'  would  be  a  linear  function  of  the  mechanical  load 
or  crack-opening  displacement.  Thus, 
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(9) 


Atlp  "  Atip  8ln<wt  +  ♦) 

where  w  is  the  frequency  of  the  loud  variation,  and  $  is  the  phase  delay 
between  load  variation  and  film  rupture* 

Equation  (9)  suggests  a  second  method  of  evaluating  the 
electrochemistry  of  the  crack  tip  environment,  using  the  transmission 
line  model*  Under  a  steady  electrical  bias  (a  small  DC  voltage  more 
positive  than  the  free  corrosion  potential),  a  DC  current  can  be  made  to 
flow  from  the  specimen  to  a  suitable  counterelectrode*  This  current 
originates  with  the  anodic  dissolution  at  the  crack  tip  and  is  given  by 

I(u)t)  -  h/R.  (10) 


l  • 

U 

! 

S  ■ 

i 


N 


where  h  -  vappne(j  “  vfree  corrosion* 

Substituting  equation  (9)  into  equation  (8)  and  then  substituting 
the  resulting  expression  for  into  equation  (10)  yields 

t|nFA°  sln(u>t  +  4>)i 

I(wt)  -  - - °  (10a) 

This  expression  for  I(u>t)  describes  the  crack  tip  current  expected  for 
cyclic  load  conditions* 

This  can  be  viewed  as  a  sinusoidal  current  source  originating  at 
the  crack  tip.  The  current  source  is  situated  where  the  Interfacial 
reaction  resistance  (Rr  in  Figure  2)  is  placed  in  Ztlp  In  Figure  1. 

Alternatively,  we  csn  define  an  electrochemical/mechanical 
Impedance, 

Zem  "  Z'1  (H> 

which  may  be  examined  as  a  function  of  the  frequency  (and  amplitude)  of 
the  mechanical  load  variation* 
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The  advantage  of  an  electrochemical/mechanical  impedance  study  ia 
that  the  perturbation  is  applied  exactly  at  the  point  of  interest  (the 
crack  tip).  In  the  alternative  method,  a  variable  frequency  electrical 
potential  variation  is  applied  to  the  mouth  of  a  crack  in  a  specimen 
cycled  with  a  constant  frequency  mechanical  load;  however,  this  results 
in  a  perturbation  that  decreases  in  amplitude  at  the  crack  tip  as  the 
crack  increases  in  length. 
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Section  3 


EXPERIMENTAL  METHODS  AND  RESULTS 

3.1  TEST  MATERIALS  AND  ENVIRONMENT 

Two  test  materials  were  used  In  this  program:  a  titanium  alloy  (Tl- 
6A1-4V)  and  a  medium-strength  low-alloy  steel  (HY  80).  These  materials 
were  selected  because  of  their  possible  apollcatlon  In  submarine  hull 
structures  and  for  the  range  of  ele^'rochemlcal  and  corrosion-fatigue 
behavior  these  materials  cover.  The  chemistry  and  mechanical  properties 
of  the  T1-6A1-4V  and  HY  80  are  given  In  Tables  1  and  2,  and  the 
microstructures  are  shown  In  Figure  3.  Note  that  the  HY  80  has  a 
balnltlc  microstructure  typical  of  a  quench  and  temper  heat  treatment, 
whereas  the  T1-6A1-4V  exhibits  an  alpha-beta  structure  typical  of  the 
mill-produced  hot-rolled  and  annealed  condition. 

The  test  environment  selected  for  this  program  was  the  standard 
3.5X  NaCl  solution.  This  environment  was  created  by  dissolving  the 
appropriate  amount  of  reagent  grade  NaCl  In  deionized  water.  This 
solution  was  aerated  and  continuously  refreshed  Into  the  test  chamber  at 
a  rate  of  about  1  liter/hr. 

3.2  ELECTROCHEMICAL/MECHANICAL  EXPERIMENTS 

To  generate  the  electrochemical/mechanical  data  needed  for  the 
transmission  line  model,  we  performed  three  types  of  AC  Impedance-based 
experiments  In  this  program:  (1)  planar  (uncracked)  electrode 
experiments  to  characterize  the  bulk  electrochemical  properties  of  the 
materials  without  any  mechanical  loading,  (2)  cracked  electrode 
experiments  to  characterize  the  voltage/current  response  during  active 
loading  of  a  crack,  and  (3)  cracked  electrode  electromechanical 
Impedance  experiments  to  characterize  the  voltage  displacement  response 
during  active  loading  of  a  crack. 
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Table  1 


PROPERTIES  OF  HY80  STEEL 

(5.08-cm-thick  plate,  M11-S-16216H,  ASTM  Grain  Size  8) 


Elemen tal 

Analysis 

Heat  Treatment 

Mechanical 

Properties; 

Element 

vtZ 

Parameter 

Value 

C 

0.17 

Austentized  at 

°vs 

87.8  ksi 

Mn 

0.34 

1010aC  for  120 

°lta 

106  ksi 

P 

0.010 

minutes  and  water 

ZEtS 

23. 7Z 

S 

0.10 

quenched; 

ZRA 

66. 7Z 

Si 

0,23 

Tempered  at 

CVN 

95  ft-lb 

Cr 

1.56 

682*C  for  120 

at  -120°F 

Ni 

2.88 

minutes  and  air 

Mo 

0.24 

cooled. 

Cu 

0.17 

Ti 

0.001 

V 

0.001 

definitions: 

<Jys  -  0.2Z  offset  yield  strength. 

Quta  *  ultimate  tensile  strength. 

ZE  -  percent  total  tensile  elongation  to  failure. 

Z  RA  ■  reduction  in  area  of  necked  region  at  fracture. 

CVN  »  charpy  V-notch  impact  energy;  CVN  was  determined 

under  the  following  conditions:  Type  A,  notch  V,  full  size, 
longitudinal  orientations. 
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Table  2 


The  methods  used  and  the  results  obtained  for  each  of  these  three 
experiments  are  discussed  in  the  following  subsections*  However,  before 
discussing  these  experiments,  It  Is  Important  first  to  review  some  of 
the  advantages  and  limitations  encountered  when  applying  these 
electrochemical  techniques. 

To  calculate  accurately  the  many  electrochemical  parameters  that 
can  be  obtained  from  AC  measurements,  It  Is  necessary  to  obtain  high- 
precision  Impedance  measurements  over  a  wide  range  of  frequencies, 
particularly  at  very  low  frequencies.  The  unavailability  of  adequate 
Instrumentation  has  been  the  principal  limitation  to  the 
characterization  of  cracked  or  porous  systems  by  AC  techniques  and  is 
largely  responsible  for  the  only  modest  success  of  Broers  and 
coworkers. More  recently,  Mund  et  al.^  analyzed  the  performance 
of  hydrogen  electrodes  In  alkaline  fuel  cells  systems  from  AC  impedance 
measurements  applied  to  a  transmission  line  model. 

The  current  work  of  Hills  and  coworkers*2-^  provides  good 
verification  for  the  application  of  the  transmission  line  model  and  the 
extension  of  equilibrium  AC  measurements  for  predicting  the  steady-state 
electrochemical  parameters.  In  one  set  of  experiments,  the  interfacial 
Impedance  components  were  calculated  from  the  measured  impedance  of  a 
packed  bed  of  graphite  spheres  using  the  transmission  line  model,  and 
the  results  were  compared  with  the  lnterfaclal  Impedance  measured 
directly  for  a  single  sphere. 

Figure  4  shows  the  Impedance  data  obtained  for  a  packed-bed  elec¬ 
trode  superimposed  on  the  least-squares  beat  fit  of  data  to  equation  (4). 
The  tabulated  values  of  the  double-layer  capacitance,  Warburg  coef¬ 
ficient,  and  exchange  current  density  for  the  packed  bed  were  calculated 
from  this  regression  fit.  These  values  were  compared  with  the 
lnterfaclal  Impedance  parameters  measured  directly  for  a  single  graphite 
sphere  Isolated  from  the  bed.  The  substantial  agreement  between  the  two 
sets  of  calculated  Impedance  parameters  Is  strong  verification  for  the 
applicability  of  the  transmission  line  model  to  characterizing  porous  or 
cracked  electrodes. 
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FIGURE  4  MEASURED  AND  CALCULATED  IMPEDANCE  SPECTRA 
FDR  A  BED  OF  600-^m  GRAPHITE  SPHERES 

(From  Reference  4) 


This  previous  research  clearly  shows  that  many  Important  operating 
parametero  of  porous  and  cracked  electrodes  are  contained  within  the 
measured  Impedance  dispersion.  However,  two  factors  are  fundamental  in 
the  application  of  this  method  to  the  in  situ  determination  of  electrode 
characteristics.  First,  impedance  measurements  must  be  of  sufficient 
precision  (better  than  0.1Z  for  each  impedance  component),  and  second, 
the  measurements  must  be  made  over  a  frequency  range  wide  enough  (10~^ 
to  104  Hz)  to  adequately  resolve  the  impedance  spectrum  and  to  allow 
subsequent  deconvolution  to  obtain  the  electrochemical  parameters. 

The  requirements  for  a  high-speed,  high-precision,  impedance 
measurement  system  capable  of  operating  at  very  low  frequencies  have 
been  presented  in  detail  by  Macdonald  and  McKubre.  *  In  general,  these 
requirements  are  in  conflict.  Maximum  acquisition  speed  is  available 
from  noise  or  transform  methods,  but  these  methods  are  not  very 
precise.  High  precision  is  best  achieved  with  a  frcquency-by- 
frequency  measurement  technique  such  as  phase-sensitive  detection,  but 
the  measurements  are  time  consuming.  Fortunately,  this  is  not  a  serious 
hindrance  when  such  equipment  is  automated  and  the  experiments  are 
computer-controlled.  We  used  such  automated,  computer-controlled 
equipment  in  the  experiments  described  in  the  following  subsections. 

3.2.1  Planar  Electrode  Experiments 

Experlmenta  were  performed  on  planar  uncracked  electrodes  of  HY80 
and  T1-6A1-4V  exposed  to  aerated  3.5  wtX  NaCl  to  obtain  the 
electrochemical  parameters  needed  for  developing  the  DC,  crack 
propagation  model  and  to  aid  in  interpretating  the  resulta  of 
electromechanical  impedance  testing.  The  experimental  setup  for  these 
teats  la  shown  in  Figure  5.  Specimens  were  mounted  with  one  planar 
surface  exposed  to  the  electrolyte. 

Figure  6  presents  the  Tafel  plot  obtained  for  TI-6A1-4V  at  very  low 
sweep  rates;  the  aweep  was  initiated  at  the  cathodic  limit  of  -500  mV 
versua  the  saturated  calomel  electrode  (SCE).  The  apparently  well- 
behaved  Tafel  law  behavior  evidenced  by  the  region  of  linear 
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FIGURE  5 


EXPERIMENTAL  CONFIGURATION  FOR  AC  AND  DC  TEST 
ON  PLANE  PARALLEL  ELECTRODES 


I IG&2S 


LOG  CURRENT  DENSITY  {A  crrf^) 


JA-4333-12 

FIGURE  6  TAFEL  PLOT  FOR  A  PLANE  PARALLEL  E LECTROOE 
OF  Ti-6AI-4V  EXPOSEO  TO  3.5  wt%  NaCI  ANO  25’’C, 
VOLTAGE  SWEEP  RATE  »  0.01  mV  s  1 
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relationship  between  voltage  and  the  logarithm  of  current  auggests 
charge-transfer  kinetic  rate  control  In  both  the  anodic  and  cathodic 
regions. 

The  AC  Impedance  behavior  of  this  electrode,  however,  strongly 

Indicates  that  the  reaction  kinetics  of  the  T1-6A1-4V  electrode  surface 

Is  dominated  by  the  presence  of  a  semiconducting  oxide  film  such  as 

T102*  Evidence  of  the  dominance  of  semiconductor-electrochemical 

effects  can  be  seen  in  the  Mott-Schottky  plot  presented  as  Figure  7. 

The  reciprocal  of  the  measured  capacitance  squared  Is  plotted  versus 

potential  at  25-mV  steps  from  -200  to  -1200  mV  versus  SCE.  Regions  of 

linearity  in  such  a  plot  Indicate  the  presence  of  a  space-charge 

layer  (band  bending  within  the  oxide  film).  The  type  of  frequency 

dispersion  ahown  in  Figure  7  (that  is,  the  lines  at  each  frequency  do 

not  superimpose)  is  unusual  and  may  be  due  to  impedance  effects  within 

the  electrolyte  phase  (the  electrical  double  layer)  or  to  heterogeneity 

within  the  oxide  itself.  The  extrapolated  intercept  of  this  family  of 

2 

lines  with  the  axis  at  C  +  ®  (1/C  •+■  0)  Indicates  the  position  of 
minimum  band  bending  (the  "flat-band"  potential),  which  Is  a 
characteristic  of  each  semiconducting  oxide. ^hls  intercept 
of  ~  1100  mV  SCE  Is  within  100  mV  of  that  expected  for  pure  TIO2.  This 
100  mV  offset  Is  a  result  of  the  defect  structure  and  doping  of  the  Ti- 
6A1-4V  film  that  forms  In  the  chloride  environment. 

The  HY80  specimens  showed  a  more  conventional  bulk  electrochemical 
response  than  the  T1-6A1-4V.  The  HY80  specimen  did  not  show  Mott- 
Schottky  behavior,  which  is  not  surprising  because  the  corrosion  film 
that  develops  Is  poroua  and  without  uniform  composition  through  its 
thickness;  that  la,  the  donor  density  varies  through  the  film 
thickness.  Figure  8  shows  the  Impedance  locus  or  Nyqulst^  response  of 
a  plane  parallel  electrode  of  HY80  exposed  to  aerated  3.5  wtX  NaCl,  at 
the  free  corrosion  potential.  The  depressed  semicircular  form  Is  that 
expected  for  a  metal  electrode  actively  corroding  under  mixed 
kinetic/mass  transport  control.  Figure  9  presents  data  collected  at  the 
same  electrode,  covering  a  range  of  potentials  more  cathodic  than  the 
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FIGURE  7  MOTT-SCHOTTKY  PLOT  FOR  A  PLANE  PARALLEL  ELECTROOE 
OF  Ti-6A!-4V  EXPOSED  TO  3.5  wt%  NaCI  AT  25*C 
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FIGURE  8  IMPEDANCE  LOCUS  PLOT  FOR  A  PLANE  PARALLEL  ELECTRODE 
OF  HY80  EXPOSED  TO  3.5  wt%  NaCI  AT  25’C 


IMAGINARY  COMPONENT,  B  (ms) 


free  corrosion  potential  of  HY80  (-  -600  mV  versus  SCE  In  aerated  3.5 
wt2  NaCl).  These  data  are  plotted  as  the  reduced  admittance,  B,  versus 
the  square  root  of  the  angular  frequency,  u>^^, 


B  -  -X/[(R  -  R  )2  +  X21 

3 


(12) 


1/2  1/2 

where  u>  ■  (2nf)  ,  X  Is  the  imaginsry  component  of  the  Impedance, 

R  Is  the  real  component  of  the  Impedance,  Rg  Is  the  uncompensated  series 
electrolyte  resistance,  and  f  Is  the  frequency  of  the  Impedance 
measurement  In  Hz. 

The  family  of  straight  llnea  In  Figure  9  Indicates  thst,  at 
potentials  more  negative  than  (cathodic  to)  the  free  corrosion 
potential,  the  Impedance  of  an  HY80  surface  Is  dominated  by  dlffuslonal 
effects  within  the  electrolyte.^2  At  potentials  higher  than  -700  mV, 
however,  deviation  from  a  linear  B  versus  w*/2  dependence  occurs,  as 
shown  by  the  data  at  -550,  -450,  and  -400  mV,  plotted  In  Figure  10. 

These  data  are  replotted  In  Figure  11  versus  u>.  The  straight  lines 
shown  are  regression  fits  to  the  data  points;  this  linearity  Indicates 
that  the  Imaginary  part  of  the  admittance  Is  dominated  by  a  capacitance 

a  /  •j 

of  approximately  20  pF  cm*.4  The  slope  of  dB/du)  is  equal  to  the 
double-layer  capacitance  and  can  be  seen  In  this  potential  range  because 
of  the  absence  of  dlffuslonal  control  effects. 

At  -300  mV  versus  SCE,  the  electrode  admittance  Increases  greatly 

1  /2 

snd  again  displays  llnesrlty  versua  u>i7  ,  as  shown  In  Figure  10.  Thia 
is  the  tranapaaslve  region  for  HY80,  and  the  electrode  Impedance  becomes 
dominated  by  the  maaa  transport  of  corrosion  products  from  the  electrode 
surface . 


In  summary,  the  recults  of  planar  electrode  experiments  on  bulk 
T1-6A1-4V  and  HY00  specimens  ahowed  that,  over  an  extended  potential 
range,  the  T1-6A1-4V  specimen  behsved  ss  If  It  were  covered  by  a  dense 

2 

semiconducting  oxide  layer.  It  exhibited  the  characteristic  linear  l/C 
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FIGURE  10  REDUCED  ADMITTANCE  VERSUS  cj1/?  FOR  A  PLANF  PARALLEL 
ELECTRODE  OF  HY80  EXPOSED  TO  3  5  wt%  NaCl  AT  25*C. 
ANODIC  POTENTIAL  RANGE 


A  % 


versus  voltage  behavior.  On  the  other  hand,  the  dependence  of  the 
Impedance  response  on  voltage  for  HY80  was  more  complex  although  more 
typical  of  a  corroding  electrode.  At  potentials  less  than  and  equal  to 
the  free  corrosion  potential,  the  electrode  Impedance  response  was 
dominated  by  diffusion,  presumably  of  reactant  (dissolved  O2)  to  the 
electrode.  At  potentials  up  to  200  mV  more  positive  than  the  free 
corrosion  potential.  In  the  potential  range  where  a  pseudo-passive  film 
Is  expected  to  form,  the  Impedance  appears  to  be  dominated  by  an 
essentially  voltage-independent  capacitance,  almost  certainly  due  to  the 
electrical  double  layer  In  the  electrolyte.  At  still  more  positive 
potentials.  In  the  transpassive  range,  the  electrode  Impedance  appears 
to  be  dominated  by  the  diffusion  of  reaction  products  from  the  electrode 
surface. 


3.2.2  Mechanically  Loaded  Cracked  Electrode  Experiments 

Ve  developed  an  experimental  facility  designed  to  obtain  data  on 
crack  electrochemistry  during  CF  and  SCC  of  various  materials  In  aqueous 
solutions  at  room  temperature  and  pressure.  Figure  12  diagrams  the 
experiment  arranged  for  Impedance  measurement  of  specimens  under 
mechanical  load  conditions.  Photographs  of  the  setup  are  shown  in 
Figure  13  (electronics).  Figure  14  (test  chamber),  and  Figure  15  (test 
specimen) . 


As  detailed  In  Figure  12,  Impedance  measurements  were  obtained 
using  a  Solartron  Model  1250  frequency  response  analyser  (FRA)  operated 
under  microcomputer  control.  The  Model  1250  Is  an  automated,  digitally 
demodulated,  stepped  frequency,  impedance  meter  capable  of  0.01% 


precision  for  each  Impedance  component  In  the  frequency  range  10”^  to 


6  x  10*  Hr.  This  Instrument  and  the  preceding  series  of  frequency 


response  analyzers  (Solartron  Model  1172  and  1174)  have  been  used 


extensively  In  recent  years  to  measure  impedances  In  electrochemical 
48-52 


systems . 


The  FRA  was  operated  under  the  control  of  an  Apple  II-Plus  micro¬ 
computer  equipped  with  64  kbyte  of  Internal  random  access  memory  (RAM), 
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FIGURE  12  SCHEMATIC  OF  IMPEDANCF  MEASUREMENT  SYSTEM  FOR  STRESS  CORROSION 
CRACK  AND  CORROSION  FATIGUE  TESTS 
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F ACM  ITY  FOR  MEASURING  IMPEDANCE  SPECTRA  OF  STRESS 
CORROSION  CRACKING  OR  CORROSION  FATIGUE  TEST  SPECIMENS 


COD  Gage 


EDM  Notch 


320  kbyte  of  external  RAM  (Axlon  "Ramdisk  320"),  and  a  20-mbyte  rigid 
disk  (Corvus  Model  20  MB/M).  The  control  program  was  written  at  SRI  In 
Basic  and  6502  Machine  languages  and  Is  Included  as  Appendix  A. 

Some  experiments  were  performed  under  potentlostatlc  control,  using 
a  Princeton  Applied  Research  Model  173/276  potentlostat ;  however,  most 
of  the  experiments  were  performed  under  pseudo-galvanostatlc  control. 

In  the  latter  mode,  a  resistance  several  orders  of  magnitude  larger  than 
the  cell  Impedance  Is  placed  In  series  with  the  counter  electrode,  and  a 
large  AC  voltage  Is  applied.  To  a  good  approximation,  the  current  Is 
determined  by  the  AC  voltage  divided  by  the  series  resistance. 
Independent  of  cell  Impedance. 

For  experiments  In  which  the  load  Is  cycled,  the  periods  of  the 
mechanical  and  electrical  cycles  must  be  synchroninzed  by  Integrating 
the  Impedance  data  over  one  complete  cycle  of  the  mechanical  load. 

Thus,  the  results  represent  an  average  over  a  full  cycle  of  crack¬ 
opening  displacement.  If  these  cycles  are  not  synchronized,  the  data 
may  contain  a  large  amount  of  systematic  error  and  will  Inevitably  show 
considerable  scatter.  For  example,  if  the  Impedance  measurements  are 
Integrated  over  a  time  that  Is  shorter  than  the  mechanical  load  cycle, 
multiple  determinations  (at  different  electrical  perturbation 
frequencies)  may  be  made  during  a  single  mechanical  cycle.  These 
determinations  will  cover  the  range  from  maximum  to  minimum  crack 
opening  and  will  appear  as  scatter  in  the  data  set. 

Because  It  Is  necessary  to  integrate  over  at  least  one  cycle  at  the 
electrical  perturbation,  some  error  will  be  Introduced  at  low  frequen¬ 
cies  when  the  Impedance  data  are  integrated  over  more  than  one  (but  not 
an  Integer  number  of)  mechanical  cycle.  This  error  is  Inevitable  and 
results  In  aooe  scatter  in  our  data  at  very  low  frequencies  (<  0.01  Hz). 

An  environmental  chamber  was  constructed  that  attachea  to  a  servo- 
hydraulic  mechanical  testing  machine.  This  chamber  and  the  associated 
attachments  allow  us  to  circulate  aerated  or  deaerated  3.5Z  NaCl 
solution  (or  other  electrolytes)  at  controlled  temperatures  past  our 
fatigue  precracked  compact  tension  specimens.  Tht  crack  growth  can  be 
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measured  by  optical  methods  or  by  the  compliance  method,  by  means  of  a 
water-immersible  LVDT  mounted  on  the  specimen.  A  special  composite 
reference/counter  electrode  incorporates  a  silver/silver  chloride 
electrode  in  a  0.25-inch-diameter  Lugin  probe,  positioned  about  0.1  inch 
from  the  start  of  the  fatigue  precrack.  The  fatigue  precrack  is  at  the 
bottom  of  a  2.00-inch-deep  notch  (1.50  inches  wide  by  0.010  inch  thick), 
and  the  tip  of  the  0.25-inch-diameter  lugin  probe  is  centered  at  the 
bottom  of  this  notch  to  minimize  "edge  effects”  on  the  electrochemical 
measurements.  The  counter  electrode  is  located  in  a  chamber  to  which 
the  Lugin  probe  is  attached.  This  arrangement  allows  us  to  polarize 
only  the  area  at  the  lugin  probe  tip  and  thus  measure  the  response  of 
this  area. 


The  cracked  electrochemical  AC  Impedance  measurements  were  obtained 
on  T1-6A1-4V  in  aerated  3.5%  NaCl  solutions.  These  tests  were  performed 
under  load  control  with  AK  varying  from  10-30  ksi  /in.  with  R  values  of 
0  to  0.2.  No  experiments  of  this  type  were  performed  on  HY-80. 

Figures  16  and  17  show  Impedance  data  for  the  growth  of  a  crack  in 

T1-6A1-4V.  The  data  designated  as  "no  crack"  in  Figure  16  show  the 

—3  4 

impedance  response  in  the  range  10  <  f  <  10  Hz  for  a  specimen 

subjected  to  a  5-Hz  cyclic  load,  but  showing  no  superficial  indication 
of  having  cracked.  In  the  Nyquist  domain  this  data  set  appears  as  a 
semicircle  whose  center  is  located  above  the  real  axis.  The  low- 
frequency  intercept  on  the  real  axis  is  very  close  to  2000  Q.  The  high- 
frequency  region  is  sho'm  in  an  amplified  form  in  Figure  16(b).  A 
second  semicircular  region  is  observed  with  a  high-frequency  intercept 
on  the  real  axis  at  52.8  Q.  The  high-frequency  intercept  corresponds  to 
the  series  resistance  due  to  the  electrolyte  path  between  the  tip  of  the 
reference  electrode  and  the  base  of  the  notch  in  the  specimen. 


In  the  Bode  plot  (Figure  16(c)],  these  data  exhibit  a  single 
minimum  in  the  phase  occurring  at  approximately  0.15  Hz.  The  magnitude 
information  shows  two  limiting  plateaus  (corresponding  to  the  Intercepts 
with  the  real  axis  in  the  Nyquist  plot)  separated  by  linear  regions  of 
slope  -1/2,  -1,  and  -1/4  with  increasing  frequency. 


35 


KvNsVlM 


Crack  a  =  41.3  mm 


-1000 


Increasing 

Frequency 


Increasing 

Frequency 


No  Crack  a  =  35.2  mm 


-2000 


AK  =  22  ksiVTnT 


R-0.t0Kmin/Kml< 


Load  Cycle  =  5  Hz 
Z  -  R  +  jX 


(a)  Nyquist  Plot 


JA-4333-7A 


FIGURE  16  IMPEDANCE  RESPONSE  OF  THE  Ti-6AI-4V  SPECIMEN  BEFORE 
AND  AFTER  THE  APPEARANCE  OF  A  MACROSCOPIC  CRACK 
IN  AERATED  3.5%  NaCI  SOLUTION 
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FIGURE  17  THE  EFFECT  OF  LOAD  ON  THE  IMPEDANCE  RESPONSE 
OF  A  MATURE  CRACK  IN  TI-6AMV  IN  AERATED 
3  5%  NaCI  SOLUTION 


Examining  these  linear  regions  in  order  of  decreasing  frequency,  we 
see  that  the  region  at  slope  -1/4  corresponds  to  the  high-frequency 
semicircle  in  the  Nyquist  plot,  shown  in  Figure  16(b)*  This  region  is 
present  for  unstressed  specimens  before  cyclic  loading  and  appears  to  be 
completely  unchanged  during  the  entire  growth  of  the  crack.  We 
speculate  that  this  region  is  associated  with  diffusion  in  the  plane  of 
the  machined  notch  and  perpendicular  to  the  axis  of  the  reference 
electrode.  Although  the  region  of  slope  -1/4  has  significant  extent  in 
the  plot,  the  impedances  involved  are  very  small  and  have  little 

significance  in  a  linear  (e.g.,  Nyquist)  plot  ot  the  data.  Wc  sill  not 
consider  this  zone  further  in  this  report. 

At  intermediate  frequencies,  all  data  sets  obtained  show  a  region 
with  a  slope  of  -1.  This  region  corresponds  to  a  f*1  dependence  and  is 
the  form  expected  for  the  equivalent  circuit  in  Figure  2,  at 
intermediate  and  high  frequencies,  when  diffusion  of  metal  ions  from  the 
electrode  surface  are  not  rate  limiting.  A  region  of  slope  -1  in  a  Bode 
plot  corresponds  tc  a  semicircle  with  its  center  located  on  the  real 
axis  in  the  Nyquist  representation.  The  beginnings  of  a  semicircle  of 
great  magnitude  (on  the  scale  shown)  can  be  seen  in  Figure  16(b);  the 
parameters  of  this  aemlcircle  are  determined  by  the  interfacial  reaction 
resistance,  R^.,  and  the  double-layer  capacitance,  C  (aee  Section  2). 

Impedance  data  obtained  from  unatreased  specimens  before  the 
application  of  a  cyclic  load  display  a  f“*  dependence  (Bode  slope  *  -1) 
down  to  limit ingly  low  frequencies.  Only  under  cyclic  loading 
conditions  does  a  region  of  slope  -1/2  appear,  and  this  region  becomes 
increasingly  significant  as  a  crack  grows.  Thus,  we  associate  the 
region  of  slope  -1/2  in  a  Bode  plot  with  the  presence  of  a  crack. 

Figures  16(a),  (b),  and  (c)  each  show  two  curves:  one  labeled 

"crack"  and  the  other  labeled  "no  crack"  on  the  baals  of  visual 

inspection  of  the  aides  of  the  compact  tension  specimen.  Because  both 
—1 12 

show  a  region  of  f  '  dependence,  the  curve  designated  as  "no  crack"  in 
Figure  16(c)  must  have  at  least  a  vestigial  rupture,  asaoclated  with 
cyclic  loading  of  the  oxide  film  at  the  crack  tip.  The  curve  labeled 


"crack"  was  obtained  when  a  crack  of  about  3  mm  was  observed  at  the  side 
of  the  specimen. 

A  region  of  apparent  f  A/i  dependence  (slope  -1/2)  for  the 
Impedance  within  the  crack  Is  predicted  by  equations  (4)  and  (5).  The 

v  1 

transmission  line  Impedance,  Zfcl ,  contains  terras  In  e'  ,  which  can  be 
expanded  as  a  Taylor  series: 

eYl  -  1  +  yl  +  (y  1 ) 2 /2  !  +  (y  1 ) 3 /3 !  +  .  .  .  (13) 


For  yl  <  <  1,  the  third  and  subsequent  terms  In  this  series  can  be 
neglected.  Substituting  equation  (13)  into  equatiuu  (4),  *~e  ohraln  an 
expression  of  the  form 


*  a  +  b/y  +  additional  terras 


(14) 


where  a  and  b  are  undetermined  constants  and  the  higher  order  terms  have 
decreasing  significance. 

—1  /2 

Because  y  is  a  function  of  Z^  '  [aee  equation  (5)],  the  Impedance 
of  Z^  placed  In  a  transmission  line  will  have  a  major  component  that  has 
the  square  root  of  the  frequency  power  dependence  of  the  Interfacial 
impedance  Itself.  Thus,  for  Z^  containing  parallel  resistance  (Rp)  and 
capacitive  (Cp)  elements  only. 


1/Z,  -  1/R  +  j2nfC  (15) 

1  p  p 

and  y  will  be  a  function  of  f*^;  from  equation  (14),  Z£^  will  be  a 
function  of  f~^^,  which  is  the  form  observed. 

-1  /2 

The  preceding  discussion  Is  Important  becauae  an  f  '  dependence 
is  often  attributed  to  dlffuslonal  proceases  at  a  plane-parallel 
electrode.  ^  L  By  os thematically  fitting  the  measured  Impedance  data 
to  the  form  predicted  by  equation  (4),  we  can  clearly  see  that  the 
observed  response  Is  due  to  a  capacitive  element  contained  within  a 
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transmission  line,  rather  than  to  diffusion;  the  coefficients  obtained 
do  not  correspond  to  any  known  dlffusional  process* 

The  capacitance  under  consideration  is  simply  that  of  the 
electrical  double  layer  on  the  sides  of  the  growing  crack.  A  resistance 
appears  in  parallel  with  this  capacitance  and  is  associated  with 
faradalc  processes  on  the  crack  wall.  At  present,  we  are  unable  to 
determine  whether  this  faradalc  process  is  anodic  (metal  dissolution)  or 
cathodic  (reprecipitation  of  metal  or  hydrogen  evolution).  To  resolve 
this  and  other  important  Issues,  we  constructed  a  discrete  (finite 
element)  transmission  line  model  for  the  growing  crack  based  on  the 
electromechanical  parameters  calculated  from  impedance  measurements,  as 
discussed  in  Section  4. 

The  two  data  sets  presented  in  Figure  16  are  not  greatly  different 
in  form;  the  major  changes  observed  are  in  the  characteristic  frequency 
and  the  intercept  with  the  real  axis  at  limiting  low  frequencies.  The 
characteristic  frequency  (fQ  ■  frequency  of  minimum  phase)  moves  from 
0.14  to  0.32  Hz  between  the  first  and  second  data  sets,  and  the  limiting 
low-frequency  intercept  decreases  from  3915  to  1995  Q.  This  decrease  is 
believed  to  be  an  effect  of  crack  opening.  For  the  data  set  labeled  "no 
crack,"  the  vestigial  crack,  although  short,  is  very  tight  and  the 
Impedance  is  large.  Under  the  same  load  conditions,  the  crack-opening 
displacement  Increases  as  the  crack  grows,  and  the  Impedance  becomes 
less.  In  fact,  a  plot  of  limiting  resistance  versus  crack  length  shows 
a  minimum  at  about  1-cm  crack  length,  the  resistance  then  increasing  as 
the  crack  lengthens,  indicating  that  effects  of  crack  length  becnme  more 
important  than  the  average  crack-opening  displacement. 

The  effects  of  crack-opening  displacement  were  examined  by  varying 
the  cyclic  load  conditions.  Specimens  were  cycled  at  a  mean  load,  plus 
and  minus  the  cyclic  load  perturbation.  Figure  17  shows  the  effects  of 
varying  both  the  mean  and  the  perturbation  stress  levels  for  a  mature 
crack  (■  2  cm  long)  cycled  at  0.02  Hz.  For  these  data  (and  for  all 
Impedance  dsts  presented  in  this  report),  each  impedance  determination 
represents  an  sverage  taken  over  one  complete  period  of  the  cyclic  load 


and  thus  reflects  some  weighted  average  of  crack  opening 
displacements.  Figure  17  clearly  shows  the  effect  of  crack  opening 
(proportional  to  load)  on  the  limiting  low-frequency  resistance.  From  a 
superficial  examination,  it  appears  that  the  effects  of  the  load 
perturbation  level  are  more  significant  than  those  of  mean  load. 

However,  we  have  not  attempted  to  separate  the  effects  of  mean  and 
perturbation  stress  levels  quantitatively. 

In  summary,  the  results  cf  electrochemical  experiments  on  a  cracked 
T1-6A1-4V  specimen  under  cyclic  load  conditions  in  3.5%  NaCl  suggest 
that  the  AC  impedance  appears  to  he  due  to  Farsdale  processes  at  the 
crack  tip  and  a  transmission  line  impedance  associated  with  the  crack 
walls.  We  also  observed  that,  at  least  qualitatively,  crack  length  and 
crack  opening  displacement  appear  to  affect  the  measured  Impedance. 

3.2.3  Electrochemical/Mechanical  Impedance  Experiments 

It  is  also  possible  to  define  an  impedance  in  terms  of  the 
electrochemical  response  of  a  fractured  specimen  to  a  mechanical  input 
that  causes  crack  propagation.  One  case  in  which  this  analysis  can  be 
applied  usefully  is  for  the  electrochemically  assisted  (stress 
corrosion)  cracking  of  metals  submerged  in  an  electrolyte  and  subjected 
to  a  cyclic  mechanical  load  or  stress.  With  some  simplification,  crack 
growth  due  to  cyclical  film  rupture  at  the  growing  crack  tip,  followed 
by  metal  dissolution,  can  be  considered  as  described  below. 

Figure  18  shows  schematically  the  case  of  a  crack  growing  in  to  a 
metal,  under  constant  current  control  and  cyclic  load  conditions.  The 
loading  conditions  can  be  controlled  to  achieve  a  sinusoidal  variation 
In  the  specimen  load: 


p  ■  p  +  p  sln(ut) 


(16) 
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FIGURE  18  CT  SPECIMEN  GEOMETRY  FOR  DETERMINING  THF 
ELECTROCHEMICAL  MECHANICAL  IMPEDANCE 
FOR  THE  PROPAGATION  OF  A  CRACK 


tv 


where  p  is  the  mean  specimen  load  and  p  the  amplitude  of  the 
superimposed  sinusoidal  perturbation*  To  a  first  approximation,  the 
area  exposed  at  the  crack  tip  can  be  considered  to  be  proportional  to 
the  crack  tip  opening  displacement* 


A  -  yB(CTOD) 

(17) 

CTOD  -  K2/2Eo 

ys 

(18) 

A  -  YB  K2/2Eo 

ys 

(19) 

where 

Y  *  a  proportionally  constant  determined  by  the  mechanical 

propertlea  of  the  oxide  film  and  metal  at  the  crack  tip 

B  -  specimen  thickness 

CTOD  ■  crack  tip  opening  displacement 

E  -  Elastic  modulus  of  the  specimen 

o  ■  yield  strength  of  the  specimen 
ys 

K  *  stress  Intensity  factor  «  p/B  f(a/w)33 
where 

f(a/w)  ■  w“^2(l-a/w)“3/2  (2  +  a/w)  [0.886  +  4.64  (a/w) 

-  13.32  (a/w)2  +  14.72  (a/w)3  -  5.6  (a/w)4)  (20) 

where 

p  -  applied  load 
a  -  crack  length 
w  ■  specimen  width 

54 

and  (a/w)  Is  determined  from 
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(21) 


yl.i 

Wi 


1 

V- 


(a/w)  -  1.0012  -  49. 6U  +  23.057U2  -  323. 91U3 
+  1789. 3U4  -  351 3. 2U5 


where 

U  -  [(BE [C0D]/p)1/2  +  l]-1 

COD  -  crack  opening  displacement  measured  at  a  distance 

0.345  w  from  the  load  line  towards  the  front  face  of  the 
specimen. 

The  value  for  COD  Is  determined  from 


COD  -  p  g(a/w) 


where 


0,345  1  4*  a/y 

g(a/w)  -  (1/EB) (1  +  -a7-w— )  (t  _  JiJ  [1.6396  +  11.02(a/w 


-  6.449(a/w) 


-  31.l4(a/w)3  +  47.251(a/w)A  -  18. 343(a/w) 5 ] 


Combining  equations  (22)  and  (16),  we  obtain 

COD  “  [p  +  p  sin((Jt)]  g(a/w)  (24) 

Figure  19  shows  the  approximate*  equivalent  circuit  for  the  inside 
of  the  crack,  with  the  impedance  of  the  oxlde-passlvated  walls  appearing 
In  parallel  with  the  Impedance  element  due  to  dissolution  of  exposed 
metal  at  the  crack  tip.  The  variation  in  applied  load  described  by 
equation  (16)  results  In  a  sinusoidal  perturbation  of  the  exposed  crack 
tip  area  such  that 


Because  of  the  potential  distribution  along  the  crack,  this  equivalent 
circuit  should  more  properly  be  represented  as  a  nonuniform,  finite 
transmission  line  for  the  oxide  wall  Impedance,  with  the  crack  tip  as  a 
terminating  Impedance.  Such  a  case  Is  shown  In  Figure  1. 


Rep  *  Electrolyte  Resistance  to  Crack  Walls 
Rp  *  Charge  Transfer  Resistance  at  Passive  Walls 
Cp  *  Double  Layer  Capacitance  at  Passive  Oxide 
Zw  =  Diffusional  Impedance  at  Passive  Walls 
Zp  *  Interfacial  Impedance  at  Passive  Walls 
R#a  =  Electrolyte  Resistance  to  Crack  Tip 
Rpp  *  Charge  Transfer  Resistance  at  Exposed  Metal 
C(ip  ■  Double  Layer  Capacitance  of  Exposed  Metal 
Ztjp  *  Interfacial  Impedance  at  Crack  Tip 
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FIGURE  19  SIMPLIFIED  EQUIVALENT  CIRCUIT  FOR  A  STRESS 
CORROSION  CRACK 
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Without  going  through  the  algebra,  for  the  equivalent  circuit  shown 
In  Figure  19,  we  can  obtain  an  expression  for  the  electrochemical/ 
mechanical  Impedance  (Zem)  for  this  system,  defined  as  the  ratio  of  the 
AC  voltage  that  appears  at  the  reference  electrode,  to  the  COD  due  to 
the  sinusoidal  load  variation  under  DC  galvanostatlc  conditions: 
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tip  v  ep 
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Under  optimal  conditions,  Zem,  measured  over  a  range  of 
frequencies,  can  be  used  to  deconvolve  the  equivalent  circuit  parameters 
In  a  manner  analogous  to  that  of  electrochemical  Impedance  analysis. 

The  advantage  of  this  method,  in  the  example  given.  Is  that  measurements 
can  be  made  under  dynamic  load  conditions  and  that  the  perturbation  Is 
Imposed  at  the  point  of  principal  Interest,  neat  the  crack  tip. 


In  this  experiment  we  galvanostat Ically  Imposed  a  O.S  mA  current 
between  the  counter  electrode  and  the  crack  tip.  The  potential  was 
measured  with  the  Ag/AgCl  reference  electrode  located  at  the  crack 
tip.  The  mechanical  frequency  [w  In  equation  (24)]  was  varied  from  20 
mHz  to  10  Hz. 

Figure  20  presents  the  real  versus  Imaginary  components  of  the 
electrochemtcal/mechanlcal  Impedance  response  measured  for  an  HY80  steel 
specimen  of  geometry  shown  In  Figures  15  and  18,  Immersed  In  3.5  wtX 
NaCl.  Because  of  the  form  of  equation  (24),  these  plots  are  somewhat 
more  complex  than  a  conventional  Nyqulst  plot.  Nevertheless,  these  data 
are  amenable  to  standard  methods  of  electrical  analysis.  Ue  observe 
that  the  processes  at  the  growing  crack  tip  dominate  at  low  frequencies. 
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FIGURE  20  COMPLEX  PLANE  PLOT  OF  THE  ELECTROCHEMICAL/MECHANICAL 
IMPEDANCE  FOR  THE  PROPAGATION  OF  A  CRACK  THROUGH  HY80 
STEEL  IN  3.5%  NaCI  SOLUTION  AT  25°C  UNDER  SINUSOIDAL  LOADING 
CONDITIONS 
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can  be 


and  the  properties  of  greatest  interest,  and  ^tip» 

deconvolved  at  limitingly  low  frequencies* 

Figure  21  shows  the  logarithm  of  the  electrochemical/mechanical 
admittance  (Ygm  -  l/Zem)  plotted  versus  log  frequency,  for  a  typical 
data  set,  as  a  function  of  the  number  of  mechanical  load  cycles 
experienced  by  the  specimen  during  crack  growth.  These  data  show  some 
evolution  in  the  crack  tip  parameters  (the  low-frequency,  decending 
portion  of  the  curves),  but  considerably  more  changes  appear  in  the 
crack  wall  impedance,  reflecting  the  effect  of  increased  crack  wall  area 
with  increasing  crack  length.  Also,  note  that  the  crack  grew  about  5.7 
mm  during  this  portion  of  the  test. 

To  summarize  the  experimental  results  of  the  electrochemical/ 
mechanical  impedance  data,  we  found  that  this  technique  was  sensitive  to 
crack  wall  area  at  frequencies  above  about  1  Hz  and  more  sensitive  to 
the  crack  tip  at  mechanical  frequencies  below  0.1  Hz.  To  understand  the 
details  of  these  data  and  the  electrochemical  impedance,  a  more  detailed 
electrochemical/mechanical  model  needs  to  be  constructed.  This  model 
together  with  the  impedance  data  can  then  be  used  to  measure  the 
relevant  electrochemical  and  mechanical  parameters  that  occur  in  a 
growing  SCC  or  CF  crack.  As  described  in  tne  next  section,  this 
approach  was  implemented  by  modeling  the  cracking  process  as  a  discrete 
transmittion  line. 
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Section  4 


APPLICATIONS  OF  THE  DISCRETE  TRANSMISSION  LINE  MODEL 
TO  CRACK  GROWTH  PREDICTION 

4.1  OVERVIEW 

In  developing  a  quantitative  model  intended  to  predict  the  combined 
effects  of  electrochemical  corrosion  and  mechanical  stress  on  crack 
growth,  the  first  concern  is  to  calculate  the  chemical,  electrochemical, 
and  mechanical  parameters  at  the  crack  tip,  from  the  prevailing 
conditions  external  to  the  specimen.  The  problem  of  computing  the 
mechanical  stress  and  strains  in  the  elastic/plastic  zone  surrounding 
the  crack  tip,  from  the  externally  measured  load  and  crack  opening 

IQ  IQ  CC  CQ 

displacement,  has  been  addressed  by  numerous  authors  ’  and  is 

not  the  subject  of  this  r  port. 

Instead,  we  are  concerned  with  the  manner  in  which  the  DC  poten¬ 
tial,  the  anodic  and  cathodic  currents,  the  solution  pH  and 
conductivity,  and  the  concentrations  of  metal  oxide  and  metal  ions  vary 
with  distance  and  time  within  the  interpenetrating  electrolyte  volume  of 
a  crack  growing  in  a  metal,  submerged  in  a  bulk  volume  of  electrolyte. 

In  general,  two  distinct  types  of  models  can  be  used  to  obtain  the 
required  distributions.  A  system  of  first-  and  second-order  nonlinear 
differential  equations  can  be  used  to  describe  the  processes  of  diffu¬ 
sion  and  interfacial  reaction  and  the  static  potential  distribution. 
Inevitably,  this  system  of  equations  cannot  be  solved  analytically,  and 
numerical  methods  and  mathematical  approximations  must  be  used  to  obtain 
a  solution. 

An  alternative  approach  has  been  developed  by  McKubre  for  studying 
porous  stripping  electrodes ^ ' ^2-44  ^  transmission  line 

equivalent  circuit,  described  in  Section  2,  is  considered  to  consist  of 
a  large,  but  finite  number  of  discrete  elements.  This  type  of  model 


also  requires  a  numerical  solution,  but  offers  the  advantage  of  modeling 
most  of  the  important  physical  and  chemical  processes  by  taken  into 
account  the  full  degree  of  nonlinearity  of  electrochemical  systems.  Of 
particular  significance  is  that  many  of  the  electrochemical  parameters 
necessary  to  solve  the  discrete  transmission  line  can  be  obtained  by 
making  low-frequency  AC  measurements  close  to  the  potential  of  zero 
current  and  solving  analytically  for  the  continuous  transmission  line 
impedance. 

The  basic  features  of  the  discrete  transmission  line  model,  applied 
to  a  stressed  crack,  can  be  seen  in  Figure  1.  Briefly,  we  consider  the 
resistivities  of  the  solid  (metal)  and  electrolyte  phases  to  be  inter¬ 
connected  by  the  impedance  of  the  metal/electrolyte  interface  where 
charge  carriers  change  from  electrons  to  ions.  In  the  discrete  model 
the  electrical  equivalent  circuit  is  assumed  to  consist  of  N  discrete 
zones,  each  of  length  a/N,  associated  with  the  crack  walls,  plus  an 
additional  zone  representing  the  tip.  In  general,  the  interfacial 
impedance  of  an  arbitrary  section  of  crack  wall  will  depend  on  the  local 
potential  difference  between  the  metal  and  electrolyte  phases,  the 
nature  and  concentration  of  dissolved  species  (principally  metal  and 
hydrogen  ions  and  oxygen),  and  the  presence  and  character  of  an 
interfaciai  oxide  film. 

The  selection  of  the  tip  Impedance  determines  the  nature  of  the 
electrochemically  assisted  cracking.  We  are  concerned  here  only  with 
mechanisms  of  crack  growth  by  anodic  dissolution  following  film  rup¬ 
ture.  With  appropriate  modification,  the  model  may  nevertheless  be  used 
to  test  mechanisms  of  crack  growth  involving  de-alloying  or  hydrogen 
embrittlement . 

Because  of  the  geometry  and  possible  concentration  gradients  within 
the  crack,  the  electrolyte  resistance  at  each  element  of  the  transmis¬ 
sion  line  is  not  constant  with  time  or  distance.  Cenerally,  the  metal 
phase  resistance  is  much  less  than  that  of  the  electrolyte  and  can  be 
considered  to  be  constant. 
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The  current  thst  flows  Into  the  crack  (or  the  potential  at  the 
crack  surface  near  the  crack  tip)  was  Imposed,  In  our  experiments,  by  a 
potentlostat  with  respect  to  an  external  counter  electrode.  Under  open 
circuit  conditions  this  current  (or  potential)  Is  Imposed  by  the 
reaction  of  redox  species  on  the  exposed  surface  of  the  specimen. 

4.2  ASSUMPTIONS 

We  consider  electrocheralc.\’  t'd  chemical  reactions  occurring  at  the 
crack  tip  (t)  on  the  crack  walln  w),  and  In  the  crack  electrolyte 
volume  (v).  The  Important  reactions  are  taken  to  be  as  follows: 


Reaction 


Slte(s) 

Reaction 

No. 

t,  w 

M  +  Ma+  +  ae" 

(5.1) 

V 

Ma+  +  a/2  H20  *  M0a/2  +  aH+ 

(5.2) 

w,  t 

a/2  H2  ♦  aH+  +  ae 

(5.3) 

w,  t 

M  +  a/2  H.O  J  MO  +  aH+  +  ae" 

Z  fl  /  z 

(5.4) 

All  reactions  are  considered  to  be  reversible.  Reaction  (5.1)  accounts 
for  metal  dissolution,  which  In  conjunction  with  the  chemical  reaction 
(5.2)  accounts  for  oxide  formation  by  dissolution/precipitation. 
Resctlon  (5.3)  is  simple  electrolysis,  and  reaction  (5.4)  accounts  for 
the  direct  electrochemical  reaction  between  metal  and  surface  metal 
oxide  film. 

The  following  concepts  and  assumptions  have  been  applied  to 
simplify  the  model. 

•  Only  metal  dissolution  (reaction  5.1)  and  the  crack  tip 
Is  Influenced  by  the  crack  opening  displacement. 

•  There  Is  no  oxygen  within  the  crack  volume  and  thus 
reaction  (5.3)  Is  the  only  source  of  net  cathodic 
current . 


•  The  crack  is  considered  to  be  a  wedge  and  the  area  of 
metal  exposed  at  the  crack  tip  is  determined  by  the 
tangent  of  the  crack  opening  angle,  0* 


Tan(Q)  -  COD/a  (28) 


where  COD  ■  crack  opening  displacement,  a  ■  crack  length. 

•  In  the  electrochemical  reactions  all  symmetry  factors  (a) 
are  considered  to  be  equal  to  1/2. 

•  The  chemical  reaction  (5.2)  rapidly  achieves  equilibrium. 

•  At  the  steady  state  the  activities  of  1^0  and  H2  are 
invariant. 

•  The  activity  of  the  metal  oxide  phase  is  equal  to  the 
monolayer  surface  coverage,  0,  with  a  maximum  of  1. 

•  At  all  points  the  crack  is  so  thin  that  there  exists  no 
limitation  due  to  diffusion  perpendicular  to  the  walls. 

•  The  presence  of  an  (NaCl)  ensures  that  the  electrolyte 

phase  resistivity  does  not  change  significantly  with  metal  or 
hydrogen  ion  concentrations* 


4.3  EQUATIONS 

Anodic  dissolution  at  the  crack  tip  is  considered  to  proceed  via  the 
Butler-Volmer  equation,  with  the  area  exposed,  At,  being  proportional  to 

D. 


At  -  yCOD/a  (29) 

where  y  is  a  constant,  and 

COD  -  5  +  | d |  sin(wt)  (30) 
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These  equations  (29  and  30)  were  used  as  a  simplified  version  of 
Equations  (19)  and  (24)  to  reduce  the  model  run  time*  Equations  (19)  and 
(24)  will  be  Incorporated  In  the  future.  Equation  (30)  considers  the 
effect  of  sinusoidal  COD  loading  of  the  compact  tension  on  the  dynamics 
of  crack  growth,  and  a  model  Is  being  developed  to  calculate  the 
electromechanical  Impedance  under  dynamic  conditions  that  also  Includes 
modeling  of  the  y  factor  in  equation  (29).  In  this  section  we  consider 
only  steady-state  conditions,  and  we  define  an  effective  opening 
displacement,  D,  which  Is  a  function  of  COD,  d  and  the  sinusoidal  loading 
ra  te. 

The  current  associated  with  metal  dissolution  at  the  crack  tip,  due 
to  an  applied  overvoltage,  under  steady-state  conditions  can  be  expressed 
as 


IH/Ma+  "  ‘m/M  3+  (lrD/a)  taxPl  aRT  <4VN  '  VM/Ma+  ~  Vcp' 

-  «■>  I  m  <"n  -  VM/M  a+  -  V”  <31> 

where  AV^  Is  the  interfaclal  overvoltage  In  element  N,  F  Is  Faraday's 
constant,  T  is  temperature  In  Kelvin,  V°  Is  the  equilibrium  potential,  1° 
Is  the  exchange  current  density,  subscript  M/Ma+  refers  to  the  metal 
dissolution  reaction  (5.1),  and  V£p  Is  the  local  concentration 
overvoltage  determined  by 

V  (N)  -  -  (nF/RT)  ln[H+]  (32) 

cp 

The  standard  state  for  1°  and  V  _  Is  taken  to  be  1  molar. 

cp 

Crack  growth  Is  determined  at  discrete,  equal  time  steps.  Within 
each  time  interval  reactions  (5.1),  (5.3),  and  (5.4)  are  evaluated  on  the 
crack  walls  In  each  segment  of  the  transmission  line,  according  to  the 
local  conditions  of  concentration  and  potential. 

At  the  end  of  each  time  Interval,  reaction  (5.2)  In  each  segment  is 
considered  to  be  at  an  equilibrium  determined  by 
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[M°a/2]tH+]a 


[Ma+] [H_0]a/2 


(33) 


where  the  activity  of  water,  [H20],  is  considered  to  be  unity,  and  the 
activity  of  the  oxide  film,  [M0a/2],  is  taken  to  be  equal  to  the  frac¬ 
tional  area  coverage  with  a  maximum  of  1>0  (at  one  monolayer  coverage)* 


The  extent  of  crack  growth  in  each  time  interval,  t,  is  calculated 
from  the  volume  of  metal  dissolved  at  the  crack  tip  in  a  rectangular 
parallellpiped  whose  width  is  the  width  (B)  of  the  crack  opening  angle 
and  v  and  whose  depth  is  the  increased  crack  length,  Aa.  Thus, 


(Ca+)  <<> 


AL 


nF  MM  p  Y  0  B 


(34) 


where  hW  and  p"  are  the  average  molecular  weight  and  density  of  the 
dissolved  metal,  respectively,  and  the  other  symbols  have  been  previously 
defined. 


The  concentrations  within  each  segment  of  the  transmission  line  are 
modified  by  the  combined  effects  of  (1)  chemical  and  electrochemical 
reactions,  (2)  the  diluting  effect  as  the  volume  is  Increased  with  crack 
growth  (in  flow  of  bulk  electrolyte  from  outsldt  the  crack),  and  (3) 
diffusion.  The  effects  of  diffusion  along  the  crack  length  are  accounted 
for  at  the  same  discrete  time  intervals  as  crack  growth.  A  comprehensive 
model  for  transverse  diffusion  is  presented  in  Appendix  B.  Briefly,  we 
consider  the  propagation  of  concentration  profiles  between  nearest 
neighbors  that  initially  have  uniform,  but  different  concentrations. 
Adjacent  zones  have  different  cross-sectional  areas  and  different  volumes 
in  the  wedge-shaped  crack  profile.  Different  approximations  are  uaed  for 
long  and  short  times. 


The  final  consideration  is  of  the  thickness  of  the  “passivating’ 
metal  oxide  film.  This  film  has  an  important  role  in  localizing  the 
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metal  dissolution  process  at  the  crack  tip*  We  consider  the  metal  oxide 
to  conduct  both  electrons  and  ions*  The  electronic  conductive  properties 
are  accounted  for  by  the  solid-state  reaction  between  metal  and  metal 
oxide  [reaction  (5*4)].  The  primary  role  we  envisage  for  the  metal 
oxide,  however,  is  as  a  diffusion  barrier*  At  greater  than  monolayer 
coverages,  we  consider  that  the  direct  metal  dissolution  reaction  on  the 
crack  walls  is  impeded  by  the  diffusion  of  metal  ions  across  the  metal 
oxide  membrane.  The  passive  current  density,  ip,  is  therefore  given  by 


(35) 


where  D  is  the  diffusion  coefficient,  Cg  is  the  saturation  solubility  of 
metal  ions  in  the  metal  oxide  phase,*  Cjq  is  the  concentration  of  metal 
ions  at  the  crack  electrolyte/metal  oxide  interface,  and  y^  is  the  oxide 
film  thickness.  We  simplified  the  calculation  by  assuming  that  metal 
salt  will  precipitate  when  the  activity  of  metal  ions  reaches  unity 
within  the  oxide  phase.  We  can  therefore  define  the  maximum  passive 
current  density  (for  Cg  -  1  and  Cjj  -  0)  as 

i°  -  -nro/yN 

and  the  passive  dissolution  current  in  each  segment  is  given  by 

JN  ‘  ‘p  *  V'h 


(36) 


(37) 


is 

r. 


V 

H 

i 


Note  that  we  have  not  specified  whether  or  not  the  film  is 

porous.  In  the  event  that  It  is,  the  values  of  Cg  and  D  would  change, 

but  the  principle  of  a  diffusion  barrier  remains  the  same. 
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4.4  IMPLEMENTATION 


The  model  described  In  this  section  has  been  Implemented  as  a 
computer  program  written  In  BASIC  for  the  Apple  II  microcomputer.  A 
listing  Is  provided  as  Appendix  C,  and  a  schematic  flow  chart  is  given  in 
Figure  22.  Briefly,  the  program  allows  for  the  Initialization  of  the 
physical  and  electrochemical  parameters  used  to  define  crack-growth  con¬ 
ditions.  These  parameters  are  stored  In  the  array  ZM(26)  and  have  the 
meanings  shown  In  Table  3. 

The  transmission  line  Is  considered  to  consist  of  N  segments,  where 
N  <  200.  Data  are  saved  In  the  array  ZN(9,200)  for  the  metal  and 
electrolyte  phase  voltages  and  resistances,  the  Interfacial  currents  for 
the  electrochemical  reactions  M  *•  Ma+,  H2  +  2  H+,  and  M  *■  M0a/2i  and  for 
the  hydrogen  and  metal  Ion  concentrations  In  the  crack  electrolyte,  and 
the  oxide  monolayer  thickness.  In  each  segment.  The  allocation  of  these 
variables  is  shown  In  Table  4. 

Before  the  application  of  stress  or  potential,  a  vestigial  crack  Is 
considered  to  exist,  filled  with  electrolyte  at  the  bulk  pH,  with  the 
walls  covered  In  a  monolayer  of  oxide,  and  the  metal  Ion  concentration  In 
chemical  equilibrium.  The  data  array  Is  initialized  with  these  concen¬ 
trations  at  the  outset  of  computation.  The  program  then  guesses  an 
Initial  value  of  the  applied  potential  that  will  result  In  the  net  DC 
current.  The  resistances,  potentials,  and  currents  are  then  calculated 
In  each  mesh  of  the  transmission  line,  progressively,  from  the  reference 
electrode  to  the  crack  tip.  If  the  sum  of  the  calculated  interfacial 
currents  Is  not  equal  to  the  net  DC  current  that  enters  the  crack  (within 
the  specified  tolerance),  then  the  initial  value  of  applied  potential  is 
adjusted  and  the  computation  loop  reiterated. 

Following  a  successful  solution  of  the  current  and  potential 
profiles,  the  changes  in  concentration  resulting  from  chemical  and 
electrochemical  reactions,  the  crack  advancement,  and  effects  of 
diffusion  are  evaluated,  and  the  resultant  pvoflles  are  output  to  printer 
and/or  disk. 
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Table  3 


PARAMETER  ARRAY 


_ Parameter _ 

DC  current 

Specimen  width 

Specimen  height 

Specimen  length 

Crack  opening  angle 

Bulk  electrolyte  pH 

Electrolyte  resistivity 

Metal  resistivity 

Mechcnical  stress  factor,  y 

Number  of  electrons/ ion 

Computation  interval 

Number  of  elements 

Computation  tolerance 

Maximum  number  of  iterations 

Exchange  current  density, 

M/M3* 

Equilibrium  potential,  M/Na+ 

Exchange  current  density, 
M/MOa/2 

Equilibrium  potential, 

m/mo„/2 

Exchange  current  density, 

H/H+ 

Equilibrium  potential,  H/H+ 
Chemical  equilibrium  constant 
Charge  per  monolayer  of  oxide 
Charge  per  volume  of  oxide 
Metal  Ion  diffusion  coefficient 
Passive  current  density 
Mechanical  loading  frequency 


Units 

A 

cm 

cm 

cm 

radians 

Q  cm 
Q  cm 
cm 

s 

mV 


A  cm 
V 

A  cm 


V 

A  cm 

V 


C 

C 

cm 

A 

Hz 


cm 


-2 

-3 


cm 
2  s"1 


cm 


-2 


Table  4 


DATA  ARRAY 


Array  Element 
ZM(x) 

Parameter 

Units 

M  * 

0 

Solid  phase  voltage,  Vg  ^ 

V 

1 

Electrolyte  phase  voltage,  Vg  N 

V 

2 

Solid  phase  resistance,  Rg  N- 

Q  cm-1 

3 

Electrolyte  phase  resistance. 

Re,N 

Q  cm-1 

4 

Interfaclal  current,  M/Ma+ 

A 

5 

Interfacial  current,  H/H+ 

A 

6 

Interfaclal  current,  M/M0a^2 

A 

7 

Hydrogen  Ion  concentration,  (H+] 

mole  dm  ^ 

8 

Metsl  Ion  concentration,  [Ma+] 

mole  dm 

9 

Metal  oxide  thickness, 

monolayers 

This  procedure  is  repeated  at  the  interval  specified  until  the  crack 
achieves  its  desired  length. 


4.5  RESULTS 

A  series  of  computations  was  performed  in  a  preliminary  attempt  to 
verify  the  model  predictions,  using  estimated  physical  and 
electrochemical  parameters  for  HY80  steel  in  3.5  wt%  NaCl.  The 
parameters  used  are  summarized  in  Table  5.  Figures  23  through  26  show 
the  results  of  computations  that  conform  closely  to  the  initial  crack 
velocities  measured  experimentally  for  the  compact  tension  specimens.  We 
measured  the  crack  velocity  to  be  about  2.3  x  10  ^  m  s  ^  and  the  model 
predicts  an  initial  velocity  of  6  x  10-^  m  s-*  (see  Figure  21). 

For  the  first  computation  (parameters  Cl  in  Table  5),  a  crack  was 
grown  from  a  *  a0  +  0.05  cm  to  a  ■  a°  +  0.76  cm  where  a°  *  3.6  cm  at  a 
constant  crack  current  of  +1  mA.  This  computation  used  a  COD  cycling  at 
1  Hz,  which  corresponds  to  a  AK  of  ~20  ksi  /in.;  this  AK  value  did  not 
change  appreciably  as  the  crack  grew  over  this  small  distance.  Figure  23 
ahows  the  time  dependence  of  the  crack  velocity,  crack  length,  and 
potential  at  the  orifice.  The  crack  propagates  initially  at  an 
essentially  constant  rate  of  0.6  pm  s-^,  but  after  about  10  days,  the 
crack  decelerates  rapidly;  at  30  days  the  crack  velocity  is  essentially 
zero.  Under  constant  current  conditions  the  potential  at  the  orifice 
Increaaes  mono ton leal ly  with  time,  as  the  active  zone  at  the  crack  tip 
becomes  further  removed.  Because  the  crack  velocity  decreases  rapidly 
after  the  crack  has  groun  approximately  0.5  cm,  one  is  inclined  to 
attribute  this  decrease  to  a  decreasing  AK  or  correspondingly  a 
significantly  lower  crack  tip  ruptured  area  that  decreases  as  the  crack 
length  increaaes.  This  is  not  the  caae  because  the  stress  intensity  and 
area  exposed  at  the  tip  have  changed  very  gradually  and  only  by  a  small 
amount  over  the  total  range  of  crack  lengths  calculated.  It  nevertheless 
appears  that  the  model  is  predicting  a  AK  threshold  value  near  20 
ksi  / in.  This  decrease  In  crack  growth  rate  seems  to  be  due  to 

precipitation  of  oxide  along  the  crack  length. 
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Table  5 


POTENTIAL  (mV  versus  SHE) 


0.8  i  -j  0.4 


0.4  g  -j  0.2 


15  20  25  30 


TIME  OF  CRACK  GROWTH  (days) 


JA-4333-23 


FIGURE  23  DC  TRANSMISSION  LINE  MODEL  CALCULATION  OF  THE  POTENTIAL, 
CRACK  LENGTH.  AND  CRACK  VELOCITY  AS  A  FUNCTION  OF  TIME 
FOR  A  CRACK  GROWN  AT  CONSTANT  CURRENT  IN  HYBO  STEEL 


Initial  conditions  Cl  in  Table  5. 


S  V  V"  v 
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CRACK  VELOCITY  (jtm  s- 


0 


0.6 


0.8 


0.2  0.4 


m 

I 


DISTANCE  ALONG  CRACK  (cm) 


JA-4333-24 

FIGURE  24  DC  TRANSMISSION  LINE  MODEL  CALCULATION  OF  THE  POTENTIAL 
VERSUS  DISTANCE  INTO  THE  CRACK  AS  A  FUNCTION  OF  TIME  FOR 
HY80  IN  ?.b%  NaCI 

In'tlal  conditions  Cl  in  Table  5. 


I 


(a)  Interfacial  Current .  1  -*>  Ma+ 

JA-4333-25 


FIGURE  25  OC  TRANSMISSION  LINE  MOOEL  CALCULATIONS  OF  INTERFACIAL 

CURRENT.  METAL  ION  CONCENTRATION  ANO  CRACK  ELECTROLYTE 
pH  AS  A  FUNCTION  OF  OISTANCE  INTO  CRACK 
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%  OP  CRACK  LENGTH 


FIGURE  25 


(c)  Crack  Electrolyte  pH 

JA-4333-2? 

DC  TRANSMISSION  LINE  MODEL  CALCULATIONS  OF  INTERFACIAL 
CURRENT.  METAL  ION  CONCENTRATION.  AND  CRACK  ELECTROLYTE 
pH  AS  A  FUNCTION  OF  DISTANCE  INTO  CRACK  (Concluded! 


INITIAL  CRACK  VELOCITY  (/inn  s 


0.4 


Crack 

Initiation 

Potential 


POTENTIAL  (mV  versus  SHE) 

(a)  Velocity  Versus  Potential 

JA-4333-28 

FIGURE  26  DC  TRANSMISSION  LINE  MODEL  CALCULATIONS  OF  VELOCITY 
VERSUS  APPLIED  POTENTIAL.  VELOCITY  VERSUS  POTENTIAL 
(EXPANDED).  AND  LOG  (VELOCITY)  VE RSUS  APPLIED  CURRENT 
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[Ill 


till 


[iii 


POTENTIAL  (mV  versus  SHE) 

(b)  Velocity  Versus  Potential  (Expanded  Scale) 

JA-1333-29 

FIGURE  26  DC  TRANSMISSION  LINE  MODEL  CALCULATIONS  OF  VELOCITY 
VERSUS  APPLIED  POTENTIAL.  VELOCITY  VERSUS  POTENTIAL 
(EXPANDED).  AND  LOG  (VELOCITY)  VERSUS  APPLIED  CURRENT 
(Continued) 
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Figure  24  presents  profiles  of  the  potential  difference  between  the 
solid  and  electrolyte  phases  (the  potential  driving  force  for  electro¬ 
chemical  reaction)  as  a  function  of  distance  and  time.  For  this  example 
the  equilibrium  potential  for  metal  dissolution  is  taken  to  be  -200  mV 
versus  the  standard  hydrogen  electrode  (SHE);  thus,  the  baseline  of  this 
plot  represents  the  approximate  potential  of  zero  dissolution  current. 

As  expected,  the  potential  decreases  between  the  crack  orifice  and  the 
crack  tip  due  to  the  finite  conductivity  of  the  electrolyte,  the 
existence  of  electrochemical  reactions  at  the  crack  walls,  and  the 
presence  of  reaction  products  in  the  crack  electrolyte  volume.  Somewhat 
unexpectedly,  for  the  conditions  chosen,  this  potential  drop  is  not  a 
large  fraction  of  the  total  overvoltage  driving  metal  dissolution,  and 
the  potential  drop  decreases  (as  a  percentage)  with  increasing  crack 
length. 

This  observation  is  surprising  in  view  of  the  results  of  previous 

o  q 

calculations.  ‘  Nevertheless,  it  is  necessary  that  a  significant  frac¬ 
tion  of  the  externally  Imposed  potential  be  reflected  at  the  growing 
crack  tip  to  account  for  the  observation  of  a  critical  potential  for 
crack  growth. ^ 

The  origins  of  the  deceleration  in  crack  growth  rate  observed  after 
10  days  for  computation  Cl  can  be  seen  clearly  in  Figures  25(a)  through 
(c).  These  figures  show  a  normalized  crack  length,  with  the  time  of 
growth  and  actual  length  given  in  parentheses.  Figure  25(a)  presents  the 
component  of  the  total  interfacial  current  due  to  metal  dissolution  via 
reaction  (5.1),  as  a  function  of  time  and  length.  The  ordinate  is 
plotted  logarithmically  to  accommodate  the  large  range  of  interfacial 
currents.  Independent  of  time  and  length,  the  majority  of  the  dissolu¬ 
tion  current  flows  from  the  region  of  film  rupture  at  the  crack  tip. 
Despite  the  presence  of  a  passivating  oxide  film,  some  metal  dissolution 
also  occurs  on  the  walls  of  the  crack.  To  a  first  approximation,  how¬ 
ever,  the  crack  can  be  considered  as  active  at  the  tip,  partially  active 
at  the  orifice,  and  passive  (with  respect  to  metal  dissolution)  in  the 
central  segments. 
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Figure  25(b)  shows  the  concentration  of  metal  ions  (plotted  a9 
-logfmetal  ion])  within  the  crack  volume.  For  the  example  calculated  the 
concentration  is  constant  before  the  imposition  of  an  external  current, 
at  a  value  determined  by  the  chemical  equilibrium  between  metal  ions, 
hydrogen  ions,  and  the  aurface  metal  oxide.  As  the  crack  grows,  metal 
ions  become  depleted  in  the  central  aegments  due  to  the  formation  of 
metal  oxide.  Metal  ions  are  produced  by  dissolution  at  the  crack  tip, 
but  the  effect  of  this  concentration  does  not  propagate  deeply  along  the 
crack  length. 

Because  reaction  (5.2)  is  assumed  to  be  at  equilibrium,  the  deple¬ 
tion  of  metal  ions  in  the  central  segments  requires  that  the  hydrogen  ion 
concentration  also  be  low.  Thus,  Figure  25(c)  shows  a  trend  toward 
alkalinity  in  the  crack  volume,  with  local  acidification  at  the  crack 
tip. 

A  second  series  of  calculations  was  performed  to  determine  the 
effect  of  crack  current  (or  external  potential)  on  the  initial  crack 
velocity.  The  parameters  used  are  listed  as  "C2"  in  Table  5,  and  the 
results  are  summarized  in  Figures  26(a)  through  (c).  Figure  26(a)  shows 
the  steady  state  crack  velocity  for  a  vestigial  crack  0.1  cm  long, 
plotted  versus  the  potential  of  the  electrolyte  at  the  crack  orifice  for 
K  ■  10-6.  This  plot  exhibits  a  discontinuity  at  about  275  mV  versus 
SHE.  Below  this  potential  the  crack  grows  very  slowly,  whereas  at  higher 
potentials  crack  velocity  increases  roughly  linearly  with  potential.  One 
can  interpret  this  discontinuity  as  a  crack  initiation  potential  although 
a  closer  inspection  of  the  velocity  versus  potential  profile  near  200  mV 
[see  Figure  26(b)]  reveals  two  abrupt  steps  in  velocity,  one  at  180  mV 
and  the  other  at  275  mV.  Figure  26(')  shows  the  same  data  plotted  as  the 
logarithm  of  crack  velocity  versus  the  current  that  flows  into  the  crack, 
which  allows  us  to  determine  the  critical  current  for  crack  growth 
of  ~  18  aA  for  the  conditions  selected. 

The  results  obtained  for  a  chemical  rate  constant  of  10  ^  are  not 
plotted  but  shew  that,  with  larger  values  of  the  chemical  rate  constant, 
the  discontinuities  in  velocity,  shown  in  Figure  26,  disappear,  and  the 


crack  growth  rate  becomes  a  simple  monotonic  function  of  potential  (or 
current).  A  second  major  effect  of  the  chemical  rate  constant  is  that 
large  values  of  the  chemical  rate  conatant  reault  in  general 
acidification  of  the  crack  electrolyte. 


Section  5 


CONCLUSIONS 


From  the  results  of  this  program,  we  draw  the  following  conclusions 
regarding  the  use  of  AC  impedance  techniques  and  transmission  line 
modeling  to  obtain  electrochemical/mechanical  properties  states 
distributed  within  actively  growing  environmentally  assisted  cracks: 


The  AC  Impedance  response  of  the  uncracked  T1-6A1-4V  was  found 
to  behave  as  if  the  titanium  alloy  were  covered  with  a  dense 
TIOj  semiconducting  oxide  layer.  The  HY80  exhibited  more 
complex  behavior  although  more  typical  of  a  corroding 
electrode.  At  transpassive  and  cathodic  potentials  the 
electrode  Impedance  behavior  was  dominated  by  diffusional 
behavior,  whereas  in  the  potential  range  where  a  pseudo-passive 
film  is  expected  to  fom,  the  impedance  response  was  dominated 
by  the  double-layer  capacitance  In  the  electrolyte. 

The  AC  impedance  response  measured  on  cracked  T1-6A1-4V 
specimens  under  corrosion-fatigue  conditions  was  found  to  depend 
on  environmentally  induced  crack  growth,  total  crack  length,  and 
change  in  stress  intensity  (AK).  The  electromechanical 
Impedance  of  HY80  steel  undergoing  corrosion  fatigue  was  modeled 
as  a  simplified  equivalent  circuit.  The  electromechanical 
Impedance  HY80  was  also  measured.  This  electromechanical 
impedance  t  hnique  was  sensitive  to  crack  wall  area  at 
frequencies  .oove  1  Hz  and  more  sensitive  to  the  crack  tip  at 
mechanical  frequencies  below  0.1  Hz. 


A  DC  transmission  line  model  of  growing  corrosion-fatigue  cracks 
constructed  from  HY80  undergoing  corrosion  fatigue  was  used  to 
predict  crack  growth  rates  and  potential  and  chemistry 
distributions  within  the  crack.  The  predicted  crack  growth 
rates  were  In  reasonable  agreement  with  the  experimentally 
measured  values.  The  model  also  predicted  a  reasonable  AK 
threshold  for  crack  growth.  At  present,  the  model  treats  the 
crack  tip  region  rather  simply,  but  several  models  exist  that 
could  be  Incorporated  rather  easily  In  the  future. 


The  results  of  this  program  have  shown  that  a  technique  combining  AC 
Impedance,  electrochemical/mechanical  Impedance,  and  transmission  line 
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modeling  offers  the  potential  to  determine  both  the  mechanical  and 
electrochemical  processes  that  occur  within  streaa  corrosion  and 
corrosion-fatigue  cracks*  At  present,  it  is  one  of  the  few  techniques 
with  this  capability.  The  results  did  not  conclusively  show  that  the 
combination  of  AC  impedance,  electromechanical  impedance,  and 
transmission  line  modeling  can  be  uaed  quantitatively  to  predict  crack 
growth  ratea  because  the  transmission  line  model  was  not  used  to 
reproduce  the  electromechanical  impedance.  We  recommend  additional 
effort  to  use  the  transmission  line  model  in  conjunction  with  the 
electromechanical/electrochemical  impedance  to  obtain  the  critical 
mechanical  and  chemical  parameters  at  the  crack  tip.  It  is  clear  that, 
with  more  modeling  and  experimental  work  on  a  wider  variety  of  materiala, 
load  cyclea,  and  aqueous  environments,  the  electrochemical/mechanical 
technique  developed  in  this  program  can  become  an  important  research 
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Appendix  A 


IMPEDANCE  DATA  MANAGEMENT  AND  CONTROL  PROGRAM 


1  REN  *0NSS.NA1N  15/JUN/19B3 

3  REM  B V  H.C.H.MCKUBSE 

10  HINElii  36300:  LONEH;  24577 

90  ZNZ  =  120:  IF  NPI  <  >0  THEN  1000 

100  01 H  ZH<12B) , ZN (5, ZNZ) , ZZ (<7) ,ZT (32) ,Z0* (fl) 

110  PRINT  CHRI  (4);*BL0A0  B.FRA* 

140  ZN(1)  -  1: ZN (2)  =  lrZIKll)  =  100: Z«(12)  =  1:ZM(13)  =  1:ZH(21)  =  1:ZNI22)  =  1;ZM<23)  =  1:ZN(24)  =  1: 
ZNI25)  =  1 : ZH(34)  =  65535:ZH(35)  =  1:ZN(3B)  =  ZNI:ZH(39)  =  1.1:ZHI44)  =  I:ZHI45)  =  020l:ZN(49)  =  2 
00 

150  ZH (51 )  *  1 : ZN (53)  =  1:ZM(55)  =  1:ZHI56)  =  ZNI:ZI1I57>  =  1:ZN(B4)  =  2:ZH(B5)  =  2:ZMIB6)  =  2:ZN(121)  = 


2:ZH(122)  =  l:ZN(123)  =  1 
160  X  =  FRF  (0):SBI  =  6:0BI  =  1: VBZ  =  23 
170  IF  SZ  <  4  OR  SZ  >  6  THEN  SI  =  5:01  =  1:VZ  =  0 
200  60SUB  20000 

300  DATA  ■FRVANVBlVNFVMVSYVLEVSRVSLVACVASVAvyAEVAPVOAVSSVMVSflVS 
G'.'RG1 

320  OATA  ,IS“.,IC,,"nS,1,MC,.,HA*,“AU-,,RA,,,OE,.,OC"."IP,,"SA,,“Sl,,*RE,,,NA*t ,HI-t,60,(,60,f,8S-1,BR 
VNHVNSVSS  VHSVSCVSOVCO* 

340  OATA  *VE* , “AA* , ’PO" , ' 1 1", *  XH“, " XL1 , 1 XO* "VL", ' TO", "SP", ■PL* , "SM" , “ST* , “PG*, “VA*, "UP", *0 
NVFDVFOVFCVLl* 

360  OATA  "BKVILVIEVICVIFVIBVIOVIIVIOVIPVEPVJPVCP VTNVTTVOPVFN* 

3B0  OATA  ■AF-.'BF-.-RF-, -TF*. *UA* , "UB* . *UR‘ , *UT* . *U6* , *LA* . "LB* . "LR" , 'LT* , 'LG" . *SV , *PP* , *PS", 'OS-^OT 
'RR', *RH' 

400  OATA  •NRVNAVPNVERVAIVARVTS' 

1000  TEXT  :  POKE  34,0:  HONE  :  PRINT 
1100  PRINT  :  PRINT  'PROGRAn  OPTIONS  ARE:* 

1120  PRINT  :  PRINT  *  1...SET  UP  SOLARTRON" 

1140  PRINT  :  PRINT  *  2... HELP* 

1160  PRINT  :  PRINT  *  3... TAKE  HEASUREHENTS* 

11B0  PRINT  :  PRINT  *  4... CHAIN  BOOE  PLOT' 

1200  PRINT  :  PRINT  *  5... CHAIN  NYQUIST  PLOT* 

1220  PRINT  :  PRINT  *  6... SAVE  DATA  ON  DISK* 

1230  PRINT  :  PR INT  *  7.. .RETRIEVE  OATA  FROfl  DISK* 

1240  PRINT  :  PRlNr  *  B...SET  TIME* 

1245  PRINT  :  PRINT  *  9... CHAIN  ANALYSIS  PROGRAM" 

1250  PRINT  :  PRINT  *  0...EX1T  PR06RAM* 

1260  U 13)  =  0:  VTAB  (23l:ZZ(4)  =  0:ZZ(5)  *  9:  60SUB  10100 
12B0  IF  ZEZ(OI  =  1  THEN  1000 

1300  ON  ZZ (3)  GO  O  1500,2500. 3000,4000, 5000.6000.7000. B000.9700 

1320  HONE  :  VTAB  9:  PRINT  *EI1T  PR06RAH*:  PRINT  :  PRINT  ’TYPE  60T0  1000  TO  RESUME’:  END 

1500  60SUB  15000:  GOTO  1000 

1900  ZMI12B)  =  l:M*  =  *T1F2*:  60SUB  29000: ZZ<4)  *  ZZt 3):M4  *  *T1F0*:  GOSUB  29000:  PRINT  :  PRINT  *RS  =  * 


A-l 


; ZZ (4) :  PRINT  :  INPUT  “USE  TNIS  VALUE  (Y/N)  ?  *i«$:  IF  HI  <  >  “Y“  THEN  15000 
1920  ZH (51 )  -  ZZ (4) :  GOTO  1500:  REN  RP 
2000  NONE  :  PRINT  “INITIALISE1 
2020  G0SU8  20000: N*  *  “:'TT2’“:  G0SU8  21000 
2040  GOSU0  29000:  REN  KE1TNLEY 
2060  GOTO  15000 

2500  NONE  :  VTAB  5:  PRINT  '  HELP:1:  PRINT  “  — VTAB  10:  PRINT  “  TNERE  IS  NO  NELP  AVAILABLE  HER 
E.“:  PRINT  “LIFE  IS  LIKE  TNAT.  IT  CAN  SONET 1NES  8E“:  PRINT  “NARO  ANO  CRUEL.1:  PRINT  1  IF  YOU  OON 
’T  KNOW  NNAT  YOU  ARE  001NG“ 

2510  PRINT  “IT  IS  POSSIBLE  THAT  YOU  SNOULO  NOT’:  PRINT  “8E  001NG  1T.“:  PRINT  *  PLEASE  00  NOT  NES1TAT 
E  TO  NOT1:  PRINT  “CALL  NE.“:  VTAB  19:  NTAB  9:  PRINT  “HIKE1 
2520  VTAB  23:  GET  201(0):  PRINT  ZOI(O):  GOTO  1000 
3000  NONE  :  VTA8  3:  PRINT  “NANE  ANO  LOCATION  OF  SAVE  OATA  FILE' 

3010  G0SU8  11500:  PRINT  :  PRINT  “FILE  NANE  MILL  BE  PRECEOEO  8Y  DNSS(0)-“:  PRINT  :  INPUT  “ENTER  FILE  NA 
NE  :  ONSS(O)-“:i0l(O) 

3015  1  *  FRE  (0) : ZFt  =  20t(0):2ZlB)  *  77: ZOt(O)  -  “ONSS<P)-“  ♦  ZFt:  GOSUB  11400:ZF*  -  “0NSS(0)-“  ♦  ZF$ 

♦  “♦“ : ZZ (8)  =  78 

3020  PRINT  :  PRINT  “STARTING  FILE  I  “;:ZZ(3)  =  ZFX:ZZ(4)  =  0:ZZ(5)  ^  100: ZZ(6)  *  1:  GOSUB  10100:  IF  ZE 
KOI  =  1  THEN  3020 
3030  ZFT  *  ZZ(3I:ZN(3)  =•  ZN(tlB) 

3040  i  --  FRE  (0):  POKE  34,0:  NONE  :  POKE  34,3:  VTA8  1:  PRINT  “HEASURENENT  HOOE  »  *;2N(I27):  VTAB  2:  PRINT 
ZFt; ZFX:  VTA8  3:  PRINT  “SCALE  FACTOR  “;ZN(53) 

3050  Z N ( 4 7 .)  =  0:  NONE  :  VTA8  5:  PRINT  “STATUS  PAGE':  PRINT  “ . “:  VTA8  B:  PRINT  “TOTAL  I  OF  POI 

NTS  “;ZNT:  VTAB  10:  PRINT  “TINEO  SHEEP": :  HTAB  19:  PRINT  “ON*:  IF  2H(125)  *  0  TNEN  VTAB  10:  NTAB 
19:  PRINT  “OFF1 

3060  VTAB  12:  PRINT  “OC  VOLTAGE  SHEEP': :  HTAB  19:  PRINT  “0N“:  IF  ZH(I20)  :  0  TNEN  VTAB  12:  NTAB  19:  PRINT 
'OFF' 

3070  VTAB  14:  PRINT  “FREQUENCY  SHEEP*;:  HTAB  19:  PRINT  “0N“:  IF  ZN(34»  *  ZHI35)  THEN  VTAB  14:  NTA8  19 
:  PRINT  “OFF1 

3080  VTAB  16:  PRINT  “HARN0N1C  SHEEP';:  NTAB  19:  PRINT  “0N“:  IF  ZN(122)  -  ZHU23)  TNEN  VTAB  16:  NTAB  1 
9:  PRINT  “UFF“ 

3090  VTA8  18:  PRINT  “INPUT  FRONT /REAR" ; i  HTAB  19:  PRINT  “FRONT":  IF  ZN(50)  )  0  THEN  VTA8  18:  HTAB  19: 

PRINT  “REAR  “:  IF  2NI50I  >  1  THEN  VTAB  18:  HTA8  19:  PRINT  “80TH“ 

3100  VTAB  2:  HTAB  2B:  PRINT  “START  9  *;2N<124):  VTAB  3:  NTA8  29:  PRINT  “STOP  9  “; Zf1(  125) 

3110  IF  ZNU26I  >  0  THEN  G0SU8  18200:  IF  Z2<3)  ♦  ZZ (4)  /  60  <  ZNI124)  THEN  FOR  11  =  1  TO  4444:  NEXT 
11:  GOSUB  18000:  GOTO  3110 
3120  FPY  -  1:  IF  ZNi50>  -  2  THEN  FRY  =  2 
3130  ZNI56I  -  INT  (ZNY  >  <FRX  J  2N(55)  «  ZN(57)I) 

3140  ZN < 39 )  r  (ZN(34l  /  78(35)1  a  (1  /  ZNI56I) 

3145  ZNiO.O)  --  ZNt:  G0SU8  18200: ZN ( 1 . 0)  *  ZZI3)  ♦  ZZ (4)  /  60:ZN(2.0>  =  ZN U 27 ) : ZN ( 3, 0 )  *  ZN 453) t ZN < *, 0 ) 

=  1250 

3150  VTAB  23:  PRINT  2N(39I .ZN!57) .ZNI56) 

3160  VTA8  5:  CALL  -  958:  PRINT 

7170  J  =  0:  IF  ZNU2I)  =  1  THEN  ZN(3)  -  ZNUI8>:  REN  INITIAL  OC 

3175  GOSUB  3900:  REN  SET  DISPLAY 

3180  II  =  3:  GOSUB  17000 

3190  FOR  JZN  =  1  TO  441  I  ZNU17I:  NETT  JZN 

3200  2N(1)  =  ZN!35):  REN  F  NIN. 

3210  11  =  1:  GOSUB  17000 

3220  ZNI29I  *  0:  IF  ’Nil)  >  =  ZN(49I  THEN  ZN(29)  :  1:  REN  COUPLING 

3230  11  =  29:  GOSUB  1?100 

3240  ZNt 25)  =  ZNT 1 22) :  REN  INITIAL  HARNONIC 


L  /jV. 
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3250  11-25:  60SUB  17000 

3260  FRX  =  0: ZM (30)  =  0:  IF  ZHI50)  =  1  THEM  ZHI30)  =  IsFRl  ■  1:  REM  REAR 

3270  11  =  30:  60 SUB  17100: ZH(45)  =  0100:  IF  ZHI127)  =  2  THEN  ZH145)  *  0200 

32BO  IF  ZNII27I  =  0  THEN  ZM (45)  =  0201 

3290  J  =  J  ♦  1 :  IF  J  =  ZHX  AHO  ZHII2BI  =  I  THEN  3620 

3295  IF  J  >  ZN1  THEN  35B0:  REN  FULL 

3300  11  =  (5:  60SUB  17000:  60SUB  IBOOO:II  =  32:  60SUB  16900:  REN  HEASURE 
3310  HI  =  ■ : *0P2, 1  * 1 1  =  4:  60SUB  24000:  IF  ZZI4)  <  >  0  THEN  3300 
3320  FOR  ZS  1  1  TO  3:ZH(ZS  -  I,J)  =  ZZIZS):  NETT  ZS 

3330  IF  ZNI127)  >  3  THEN  II  =  45:ZHI45)  =  0200:  60SUB  17000:  60SUB  1B000:H$  =  “:'0P2, I'*: II  -  4:  GOSUB 
24000:  REN  CHI2 

3340  IF  ZH (127)  <  3  THEN  ZN(3,JI  =  1:ZH(4,J)  =  0 

3350  IF  ZHU27)  =  3  THEN  ZHI3,JI  -  ZHU  JlsZHUJI  =  ZNI2, J):ZN(I, J)  --  ZZ(2):ZN(2,J)  *  ZZ(3) 

3360  IF  ZHU27)  =  4  THEN  ZHI3.JI  -  ZZ (2)  -  ZN ( 1 . J ) : ZN (4 , J)  *  ZZ(3)  -  ZH(2,J) 

3370  IF  ZNII27)  =  3  THEN  ZHI3,J>  *  ZN(1,J) :ZN(4,JI  *  ZH <2, J) : ZN ( 1 , J)  =  ZZ 42) ; ZH42, J)  -  ZZ (3) 

33B0  IF  ZHU27)  =  5  THEN  GOSUB  3B00:  GOTO  3420 

3390  GOSUB  9500:ZN(3,J)  -  I  /  ZN(51):ZHI4,J)  =  6.2B232  )  ZH(0,J>  I  ZH(52):  GOSUB  9500:ZHU,J>  *  ZN(I,J 
)  /  ZH (53) : ZH (2, J>  =  ZN(2,J)  /  ZHI53):  REN  SCALE 
3400  GOSUB  99B0:ZHI3, J)  *  ZHI3):  REN  HA6.6PHASE 
3420  II  =  l:Ht  *  ,:"’ER,,s  GOSUB  24000:  VTAB  21:  HTAB  12:  PRINT  ZZ (II) 

3430  ZZI3)  =  7: ZZ ( 4)  =  ZH 4 122) : Z Z (5)  =  ZN(llB) :ZZ46)  =  ZHI35):  GOSUB  3990 
3432  ZZ< 3)  s  B: ZZ ( 4 )  =  ZHI123I:ZZI5I  =  ZD ( 119) : ZZ (6)  *  ZHI34):  GOSUB  3990 
3434  ZZI3)  *  9: ZZ (4)  =  1:ZZ(5)  =  ZH(120):ZZ(61  *  ZN(39):  GOSUB  3990 
3436  Z Z (3)  *  10: Z Z (4 )  -  ZH(25»:ZZ(5I  =  ZH(3):ZZ(6I  *  ZH 4 1 ) :  60SUB  3990 
3440  VTAB  12:  HTAB  9:  PRINT  J  -  IjVjZNI;  TAB1  25); J;*/';ZNl 

3450  FOR  11  =  0  TO  5:  VTAB  113  ♦  11):  HTAB  6:  CALL  -  B6B:  PRINT  ZHUI.J  -  1);  TABI  24);ZH(II,J>:  NE1T 
11 

3460  INVERSE  :  VTAB  20:  HTAB  12:Z0U0)  =  'FRONT1:  IF  ZH430)  *  1  THEN  ZOIIO)  *  1  REAR1 
3462  PRINT  ZOKO);:  HTAB  2B:Z0K0)  «  ‘0C":  IF  ZH(29)  «  1  THEN  ZOKO)  *  'AC 
3464  PRINT  ZOIIO);:  HTAB  I6:Z0I(0)  *  'OFF':  IF  ZH(14>  *  1  THEN  ZOUO)  =  'OH  1 
3466  PRINT  ZOKO):  NORNAL 

34BO  VTAB  4:  CALL  *  B6B:  IF  ZH(47)  =  999  THEN  PRINT  'TERHIHATE' 

3490  IF  ZNI50)  *  2  ANO  FRX  *  0  THEN  ZHI30)  *  1:FR1  *  1:  GOTO  3270:  REN  REAR 
3500  IF  ZH  (25)  <  Zlf  (123)  THEN  ZN125)  =  ZH(25)  ♦  1:  GOTO  3250:  REN  LOOP  H 
3510  IF  ZHU)  >  =  Zl! ( 34 )  THEN  3540 

3520  ZN(1)  -  Zlf ( 1  >  I  Zlf ( 39) :  IF  ZHU)  >  ZHI34)  THEN  ZHU)  *  ZNI34) 

3530  GOTO  3210 

3540  IF  SGN  (ZHU 20))  *  SGN  (ZHU19)  -  ZH(3)>  THEN  ZNI3I  -  ZH(3)  ♦  ZH(  120) :  6010  31B0 

3550  REH  OC  L1H1T  RE AC HE  0 

3560  IF  ZRU21)  -  1  THEN  3IB0 

3570  GOTO  3210 

35B0  REH  FULL 

3590  GOSUB  1 8000 

3600  ZH ( 1 24)  =  ZRU24I  ♦  2HU26):  REH  INC. TIRE 

361n  I  *  FRE  (0) :ZFX  *  ZFX  ♦  l-.ZOKOI  *  ZFl  ♦  STR*  (ZFI):ZZ(B)  7B:  GOSUB  11400:  REH  SAVE  OATA 
3612  IF  (ZZ (3)  ♦  ZZM)  /  1001  >  ZRU25)  OR  ZR(47)  *  999  THEN  1000:  REH  TERR1NATE 
3616  IF  ZHU2H  -  2  THEN  ZH(3)  =  ZH(3l  ♦  7HU20U  IF  SGN  IZHU2))  >■  SGN  (ZHU19)  -  ZH(3I)  THEN  1000: 
"EH  DC  UNIT 

3611  SOTO  3040:  REH  CLOSE  LOOP 

3620  ZHU)  »  60:11  *  1:  60SW  17000: ZKQ.ZHI)  «  60:  FOR  A  «  1  TO  441:  NEXT  A:M  *  'TlFl't  MSUI  29200:Z 
N(2,ZNX)  *  ZZ(3):RI  *  'TlFO';  G0SU1  29200: ZNll.ZRXI  *  ZZ(3) 


3630  ZNIIB)  -  0:11  =  18:  60SU8  17000:  FOR  II  =  1  TO  441:  HE  XT  II:HI  -  'TIF1'j  60SU8  29200: ZH(4,ZNX>  =  Z 
Z(3):H*  =  'TIFO':  60SU8  29200: ZH(5,ZNX)  -  ZZ(3) 

3640  II  =  20:  G0SU8  16900:  GOTO  3420 
3800  REH  CHI=I ,CH2=V 

3820  ZHI3.J)  =  ZH(25):ZH(4,Ji  =  ZZ(2):ZHI5,J)  -  ZZI3) 

3840  ZZI2)  =  1  /  ZH(51):Z2I3)  =  6.2B232  t  ZNIO.J)  4  ZH(52):ZZI4I  =  ZNII.J)  I  ll  12)  -  ZHI2.J)  I  ZZI3):ZH 
(2,J)  --  ZNII.J)  *  ZZ(3)  +  ZHI2.J)  I  ZZI2):ZN(I,J)  --  ZZ(4):  REH  I=E/R 
3860  ZNII.J)  =  ZHII.J)  I  ZHI53I :ZHI2, J1  -  ZNI2.J)  I  ZHI53) ;ZNI4, J)  =  ZNI4.J)  <  ZHI53):ZNI5, J)  =  ZNI5.J) 
I  ZH(53):  REH  SCALE  HULTIPLIER 
3880  RETURN 
3900  REH  DISPLAY 
3920  REH 

3930  VTA8  6:  CALL  -  958:  HTAB  B:  PRINT  *H  DCV  FREQUENCY":  PRINT  MNIT.*:  PRIHT  •FINAL":  PRINT  • 
STEP1:  PRIHT  'CURR.':  FOR  II  =  0  TO  5:  VTAB  (13  ♦  III:  PRIHT  ■fill:  HEXT  II 
3940  VTA8  20:  PRINT  'INPUT  FROH  :  C0UPLIN6  :C0HP  ' 

3960  VTA8  21:  PRIHT  'LAST  ERROR:  RETURH 

3990  VTA8  ZZ(3):  HTAB  7:  CALL  -  86B:  PRINT  ZZI 4);  TA8(  I1);ZZ(5);  TAB(  IB)jZZ(S);  RETURH 
4000  ZD*(0)  =  'BODE':  GOSUB  19900 
4020  NPX  =  1:  CALL  520"DHSS. BODE" 

5000  ZD4I0)  =  "NY0UI5T" :  GOSUB  19900 
5020  NPX  =  I:  CALL  520'DHSS.NYfiUIST* 

6000  REH  SAVE  DATA  TO  DISK 

6010  60SU8  11500:  PRIHT  :  PRINT  'FILE  NAHE  KILL  BE  PRECEDED  BY  DHSS(D)-':  PRINT  :  INPUT  'ENTER  FILE  NA 
HE  =  DHSSlDI-'iZD*lO):ZD<(0)  =  *DHSS(D1-'  ♦  ZD*(0> 

6020  ZZ 181  =  7B:  60SUB  11420:  GOTO  1000 

6500  ZDIIO)  =  ZD4I0)  ♦  STR*  (  INT  (ZZ(31  ♦  0.5)1:  RETURN 

7000  REH  RETRIEVE  DATA  FROH  DISK 

7040  I  =  FRE  (0):  GOSUB  11500:  PRINT  :  PRINT  'CURRENT  FILE  NAHE  =  ': ZF* t ' ZFX:  PRINT  :  INPUT  •  NEH 

FILE  NAHE  =  DHSS ID) ZD$ <0):ZD« (0)  *  'DHSS(D)-'  ♦  ZD*<0):  IF  ZD*(0)  (  >  'Df4SS(D)-'  THEN  ZF*  =  Z 
DI(O) 

7060  X  =  FRE  (0):  INPUT  '  HEN  FILE  *  =  *;ZFX: ZD*(0)  =  ZF*  ♦  '*'  ♦  SIR*  (ZFX) 

7080  ZZ(8)  *  78:  G0SU8  11300:  GOTO  1000 
8000  HONE  :  REH  SET  TINE 

8020  GOSUB  1B200:  VTA8  B:  PRINT  'NEN  TINE:':  HTAB  9:  INPUT  'HOURS  ';JZN:  IF  JZN  )  99  OR  JZN  <  0  THEN  8 

020 

8040  VTAB  11:  HTAB  7:  INPUT  ‘HINUTES  'jZH(83):  IF  ZHI83)  >  59  OR  ZHI83)  <  0  THEN  B040 

8060  II  *  83:  60SU8  17200:  60T0  1000 

9000  REH  DISPLAY  COUNTER  ROUTINE 

9010  HTAB  30:  PRINT  '*  JZN; '  ':  RETURN 

9240  PR$  3:  PRINT  'Jl:':  PRINT  " ZDI (0) ; '  1 PRINT  'B^JH:':  PR*  0:  INPUT  ZD*(0):  IN*  Oi  RETURN 

9500  REH  COHPLEXDIVISION 

9520  ZZ(5)  *  ZN(3.J)  *  ZNI3.J)  ♦  ZNI4.J)  I  ZNI4.J) 

9540  ZZ(6)  =  (ZN(1, J)  »  ZN13.J)  ♦  ZNI2.J)  I  ZNI4.J))  /  ZZ(5> 

9560  ZZI4)  *  IZNI2.J)  I  ZNI3.J)  •  ZNII.J)  I  ZN(4,J>)  /  ZZI5) 

9580  ZNII.J)  *  ZZ(6):ZN(2.J)  =  ZZ<4):  RETURN 
9700  ZDt(O)  -  'ANALYSIS':  G0SU8  19900 
9720  CALL  520'DHSS. ANALYSIS' 

9970  REH  HA6.  4  PHASE 

4980  ZNI3.J)  *  SDR  (ZNll.Jl  I  ZNII.J)  ♦  ZNI2.J)  I  ZNI2.JI) 

9985  ZNI4.J)  =  ATN  <ZN(2,J>  /  ZNII.J))  t  i 180  /  3.141592654) 

9990  ZNI5.J)  =  LOG  IZNI3.JI)  /  L06  (10):  RETURN 
9999  PRINT  CHRI  (131:  CHR*  I4);'SAVED«SS.HAIN':  END 
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i  A  ,*»  'J 


>  v.v.vv.v.v.v.v.vi 


10000  REN  Y/N  SU8R 

10010  ZEKO)  «  0:  PRINT  'YES  :  SET  ZD$(0):ZZ(J)  =  PEEK  (36)  -  4:  IF  ZD*<0)  =  'Y'  OR  ZD$lO)  *  CHR*  ( 

13)  THEN  ZZ(3)  «  2:  PRINT  M:  RETURN 

10020  IF  ZDKO)  =  'N'  THEN  POKE  36,ZZ(3)s  PRINT  ‘NO  ':ZZI3)  -  1:  RETURN 

10030  ZEKO)  =  I:  60SUB  10050:  RETURN 

10050  REN  ERROR  SUBR 

10060  PRINT  CHR$  (7): ZEKO)  «  1:  RETURN 

10100  REN  1/P  SUBR 

101 10  CALL  -  B6B 

10120  ZEKO)  =  0:  PRINT  ZZ (3):  INPUT  1  a;ZD$(Ol:  IF  ZDKOI  =  ■*  THEN  ZZ<7)  =  ZZ(3):  SOTO  10170 
1 0130  FOR  1  *  I  TO  LEN  (ZD$(0)):Z0»  =  HID*  (ZDKO),I,I) 

10140  IF  ZQ*  =  OR  ZB*  =  OR  ZB*  =  *  *  OR  ZB*  1  OR  ZB*  =  'E*  THEN  I0I60 

10150  IF  ZB*  <  'O'  OR  ZB*  >  '9'  THEN  10050 

I0I60  NEXT  I : ZZ < 7)  -  VAL  (ZD*(0)l 

1 0170  IF  Z Z (7)  <  ZZ (4)  OR  ZZ(7)  >  ZZ(5)  THEN  10050 

1 0180  ZZ(3)  =  ZZ <7 ) :  IF  ZZ(6)  =  I  THEN  ZZ<3)  =  1NT  (ZZ(7») 

1 0190  FOR  1  -  0  TO  200:  NEXT  1:  RETURN 

11300  REN  :  LOAD  ARRAY 

11320  POKE  9B.90:  POKE  99.ZZI8):  CALL  38152 

11330  PRINT  CHR*  (4 ) ; * BtOAD* ZD«(0> ; ■ , A* :  PEEK  (38394)  ♦  PEEK  138395)  I  256 

1 1340  RETURN 

1 1 400  REN  :  SAVE  ARRAY 

11420  POKE  98,90:  POKE  99,ZZ(B):  CALL  38104 

11430  PRINT  CHR*  (4) ; *8SAVE*ZD$(0> * , A*  PEEK  (38394)  ♦  PEEK  (38395)  I  256*,L*  PEEK  (38396)  ♦  PEEK  (3 
8397)  »  256' 

11440  RETURN 

11500  REN  ASK  SLOT, DRIVE  AND  VOLUNE 

11510  VTAB  6:  CALL  -  95B:  PRINT  'POINTERS  CURRENTLY  SET  AT  VTAB  8:  PRINT  *  SLOT  I'jSI:  PRINT  * 
DRIVE  I'sDK  PRINT  '  VOLUNE  l';VI 

11520  VTA8  12:  PRINT  'CHANSE  POINTERS  (Y/N)  ?  SET  ZD*IO):  IF  ZD*<0)  =  'N'  THEN  11550 
11530  PRINT  :  PRINT  :ZZ(3>  =  SI:ZZ(4!  *  4:ZZ(5)  =  6:ZZ(6)  *  0:  PRINT  *  SLOT  60SU8  10120:  IF  ZEI 

(0)  =  I  THEN  1 1530 
11532  SI  =  INT  <ZZ(3»  ♦  0.5) 

11535  PRINT  :ZZ<3)  =  DI:ZZ(4)  =  1;ZZ(5)  =  2:ZZ(6)  =  0:  PRINT  *  DRIVE  *;:  S0SU8  10120:  IF  ZEKO)  *  I  THEN 
11535 

11537  DI  =  INT  (ZZI3)  ♦  0.5):  IF  SI  <  >6  THEN  VI  =  0:  60T0  11550:  REN  FLOPPY 
1  1540  PRINT  :ZZ(3>  *  VI:ZZ<4)  =  1:ZZ(5)  =  67:ZZ(6)  =  0:  PRINT  •  VOLUNE  60SU8  10120:  IF  ZEKO)  *  1 
THEN  1 1540 

11542  VI  =  INT  (ZZ (3)  ♦  0.5! 

11550  PRINT  :  PRINT  'PRESS  ANY  »Ey  TO  CATAL06  ':  6ET  ZD»(0):  PRINT  :  PRINT  CHRi  (4) j'CATALOS.S'SI’.D' 


Dl'.V'VI 
11560  PRINT 


PRINT  'ANOTHER  CATAL06  Y/N  "> 

11570  IF  ZZI3I  -  2  THEN  11500 
11580  RETURN 

15000  TETT  :  POKE  34.0:  HONE  ;  VTAB  5:  PRINT  'SET-UP  OPTIONS  ARE 
15020  PRINT  :  PRINT  '  1... RETURN  TO  NAIN  NENU' 

15040  PRINT  :  PRINT  •  2. . .CREATE/NODIFY  SET-UP  FILE' 

15060  PRINT  :  PRINT  •  I.. .RETRIEVE  SET-UP  FILE  FxCN  SISK' 

15080  PRINT  :  PRINT  *  4... SAVE  SET-UP  FILE  TO  DISK* 

15100  PRINT  :  PRINT  '  5...EIECUTE  SET-UP  FILE' 

15120  PRINT  :  PRINT  '  6.. .RESET  (INITIALIZE)  SOlARTRON' 

15140  PRINt  :  PRINT  '  7. ..RESET  (INITIALIZE)  REITHLEY' 


SOSUB  10000:  IF  ZEI (01  =  1  THEN  1 1 560 
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15200  22(3)  =  I:  VTAB  (231:22(4)  =  1:22(5)  =  7:  SOSUB  10100:  IF  ZENO)  =  1  THEN  15000 
15220  ON  22(3)  SOTO  1000,15300,16600, 16640. 16700,2000, 1900 
15300  HONE  :  VTAB  5:  PRINT  ’SET-UP  HHICH  SUBSYSTEM' 

15320  PRINT  ;  PRINT  ■  I ...  GENERATOR’ 

15340  PRINT  :  PRINT  ’  2... ANALYSER’ 

15360  PRINT  :  PRINT  *  3...SNEEP(S) ’ 

153B0  PRINT  :  PRINT  ■  4...1NPE0ANCE  STANOARO’ 

15400  PRINT  :  PRINT  ’  5. .. MEASUREMENT  CONFIGURATION’ 

15460  PRINT  :  PRINT  ’  0... RETURN  TO  SET-UP  OPTIONS* 

154B0  22(3)  *  0:  VTAB  (231:22(4)  =  0:22(5)  =  5:  SOSUB  10100:  IF  ZENO)  =  1  THEN  15300 
15500  ON  22(3)  SOTO  15600,15800,16000,16300,16500 
15520  SOTO  15000 

15600  HONE  :  VTAB  3:  PRINT  ’GENERATOR’ 

15610  VTAB  6:  PRINT  ’FREBUENCV  ’;:Z2(3>  -  ZH(1):ZZ(4)  =  IE  -  5:22(5)  =  65535:72(6)  *  0:  SOSUB  10100:  IF 
ZENO)  =  1  THEN  15610 
15620  ZM ( 1 )  =  ZZ(3) 

15630  VTAB  9:  PRINT  ’AMPLITUDE  *s:ZZ(3)  *  ZN(2):Z2(4)  =  0:22(51  =  I0.23:ZZ(6)  =  0:  SOSUB  IOIOO:  IF  ZEX 
(0)  =  1  THEN  15630 
15640  ZM(2)  *  22(3) 

15650  VTAB  12:  PRINT  ’O.C.  BIAS  ’;:ZZ<3)  =  ZN(3):ZZ(4I  :  -  10.23:22(5)  =  10.23:22(6)  =  0:  SOSUB  IOIOO 
:  IF  ZEXIOI  =  1  THEN  15650 
15660  ZH<3)  -  22(3) 

15666  IF  ZM(2)  ♦  ABS  (ZM(3H  >  14  THEN  VTAB  B:  PRINT  ’CONBINEO  VOLTAGE  TOO  LMSE’:  60T0  15630 
15670  VTAB  15:  PRINT  ’WAVEFORM  (S1N=0.S8RM)  ’;:ZZ(3)  =  ZHI4):ZZ<4)  =  0:22(5)  *  l:ZZ(6)  =  I:  SOSUB  101 
00:  IF  ZENO)  =  1  THEN  15670 
156B0  2N(5)  ’  22(3) 

15690  VTAB  IB:  PRINT  ’AMPLITUDE  CONP.  (0FF=0,0NM)  ’;:22(3)  *  ZN(14):ZZ(4)  =  0:ZZ(5)  =  1:22(6)  =  I:  SOSUB 
IOIOO:  IF  ZENO)  =  1  THEN  15690 
15700  2R(I4)  =  22(3) 

15710  IF  2M ( 1 4)  *  0  THEN  15300 

15720  HONE  :  VTAB  5:  PRINT  'ANPL1TU0E  ONPRESSION  ON’ 

15730  VTAB  B:  PRINT  ’SOURCE  CHANNEL  ’::ZZ<3>  -  2Ntll):Z2(4)  =  1:Z2(5)  =  2:2Z(6)  =  1:  SOSUB  IJtOO:  IF  2 
ENO)  -  1  THEN  15730 
15740  ZN(ll)  =  1NT  (100  I  2Z(3)  ♦  0.5) 

15750  VTAB  12:  PRINT  ’VOLTAGE  LEVEL  ’{ :ZZ(3>  =  2N(I2):2Z(4)  =  IE  -  4:ZZ(5)  =  300:22(4)  *  0:  SOSUB  1010 
0:  IF  ZENO)  =  I  THEN  15750 
15760  ZN ( 1 2)  -  22(3) 

15770  VTAB  16:  PRINT  ’ERROR  PERCENT  *;:ZZ(3)  *  ZN(13):ZZ(4i  *  1:2Z(5)  *  50:22(6)  =  0:  SOSUB  IOIOO:  If 
ZENO)  *  1  THEN  15770 
157B0  ZN (13)  =  22(3) 

>5790  SOTO  15300 

15S00  HONE  :  VTAB  3:  PRINT  ‘ANAU5EP* 

158)0  VTAB  6:  PRINT  ’1NTE6RAT10N  TINE  (SECONDS)  ’;:Z2(3)  *  ZN(2)):22I4)  *  .01:22(5)  *  1E*.!Z<6)  *  0:  SOSUB 
IOIOO:  IF  ZENO)  =  1  THEN  15B10 
15920  2N<2H  *  22(3) 

15830  VTAB  9:  PRINT  ’NtASURENENT  DELAY  ’:  SOSUB  15995:22(3)  =  ZN(II):ZZ(4)  =  0:22(5)  *  IE5:ZZ(6>  *  0:  SOSUB 
IOIOO:  IF  ZENO)  =  I  THEN  15B30 
15840  ZNd!)  =  22(3! 

15850  VTAB  14:  PRINT  ’INITIAL  HARN0N1C  ’::ZZ (3)  *  2N(25):Z2'4)  »  1:22(5)  *  16:22(6)  ;  1:  SOSUB  IOIOO:  IF 
ZENO  *  I  THEN  (5850 
IjoeO  2NI25/  1  2Z(3) 
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Z5B70  VTAB  IB:  PRINT  "AUTO  INTEGRATE  OFF=0*:  PRINT  *CHI (L>=Z, CH2 (L> =2,CHZ CS) =3, CH2(S) =4": ZZ (3 )  =  ZH(2 
6) : ZZ (4)  =  0: Z Z (5)  =  4 :ZZ(6>  =  1:  60SUB  10100:  IF  ZEXIOl  *  I  THEN  15B70 
158B0  ZMI26)  -  ZZ(3> 

15900  VTAB  6:  CALL  -  95B:  VTAB  B:  PRZNT  ’INPUT  VOLTAGE  RANGE  AUTO^O’i  PRINT  “30HV=1 , 300HV=2, 3V=3, JOV 
=4.300V=5 ’ : ZZ 131  =  ZMI271:ZZ(41  -  0:ZZ(5I  =  S:ZZ(6)  =  1:  GOSUB  10100:  IF  ZEIIO)  =  1  THEN  15900 
15910  ZH127)  =  ZZ(31 

15920  VTAB  13:  PRINT  ’INPUT  COUPLING  ’:  PRINT  *0C=0, ACM. FREQUENCY  TO  CHANGE^XXX’: ZZ (3>  =  ZXX(1>:ZZ(4) 

=  0:11(5)  •-  999: ZZ 46)  =  0:  60SUB  10100:  IF  ZEIIO)  =  1  THEN  15920 
15930  ZH149)  =  ZZ131:ZH(29l  =  ZHI491:  IF  ZMI491  >  1  THEN  ZMI29)  *  0 

15950  VTAB  ZB:  PRINT  ‘INPUT  <FR0NT=0tREAR=l,B0TH=2»  ’::ZZ(3)  =  ZN(501:ZZ(4)  =  0:ZZI5)  =  2:ZZ(6)  =  1:  GOSUB 
10100:  ZF  ZEXI01  =  1  THEN  15950 

15960  Z (1(501  -  ZZ(31:ZN'301  =  ZN(50I:  IF  ZMI501  >  1  THEN  ZNI30)  -  0 
15990  GOTO  15300 

15995  PRINT  ’ (S=SECONDS. C=CVCLES>  ’ ; :  1 1  =  24:  PRINT  “CYCLES  ’;:  GET  ZDK01:  ZF  ZD4(01  =  ’S’  THEN  HTAB 
22:  VTAB  10:  PRINT  "SECONDS’;: II  =  23 

15996  RETURN 

16000  HONE  ;  VTAB  3:  PRINT  ’SHEEP (S> ... #  OF  MEASUREMENTS  *  ’;ZNX 

16010  VTAB  6:  CALL  -  95B:  PRINT  ’OC  VOLTAGE  SHEEP  IY/N1  ?  GOSUB  10000:ZH(57)  =  I:  IF  ZEXIOl  =  1  THEN 
16010 

16020  IF  ZZ<3)  =  1  THEN  ZM(UB)  --  ZM ( 3) : ZM ( 1 19 )  r  ZH(3):ZH(120)  *  0:  60T0  16100 
16030  VTAB  9:  PRINT  ’INITIAL  OC  VOLTAGE  *;:ZZ(3l  =  ZM(31:ZZ(41  *  -  10.23:ZZ(5)  *  10.23:ZZ(6)  *  0:  GOSUB 
10100:  IF  ZEXIOl  *  1  THEN  16030 
16040  ZMI11B)  =  ZZ (31 

16050  VTAB  12:  PRINT  ’FINAL  DC  VOLTAGE  *::ZZ(3I  --  ZM(119):ZZ(4)  =  -  10.23:ZZ(51  =  10.23:ZZ(61  =  0:  GOSUB 
10100:  IF  ZEXIOl  =  1  THEN  16050 
16060  ZHI1191  *  ZZI3) 

16070  VTAB  15:  PRINT  ’DC  VOLTAGE  STEP  ’::ZZ 131  =  ZM 1 1201 : Z Z (4 1  *  -  10.23:ZZ(5>  *  10.23:ZZI6)  *  0:  60SU8 
10100:  IF  ZEXIOl  *  1  THEN  16070 
160B0  ZMII20)  *  ZZ 1 31 

160B5  IF  S6N  IZMI119)  *  ZHIUBll  <  >  SGN  IZNI12011  THEM  VTAB  B:  PRINT  ’ARE  YOU  GOING  TO  CIRCLE  IMF 
1M1TY  ■>':  60T0  16030 

160BB  IF  ZMI120)  *  0  THEM  ZMI121)  *  2:ZM(571  =  1:ZHI!17I  -  0:  GOTO  16100 

16090  VTAB  IB:  PRINT  ’FILL  ARRAY  BETWEEN  DC  STEPS  IY/HI*;:  GOSUB  10000:ZMil2ll  =  ZZI3):  IF  ZHI12I1  =  2 
THEN  16100 

16093  ZMI57)  =  1MT  (1.5  ♦  IZNI119)  -  ZMIllBl)  /  ZMU201I 

16095  VTAB  20:  PRINT  ‘DELAY  BETWEEN  OCV  i/S>  ‘;:ZZ(3t  =  ZH ( 1 1 7 > Z Z ( 4 )  =  0: ZZ<5)  =  lOOO:ZZ(6l  -  0:  GOSUB 
lOlOO:  IF  ZEXIOl  --  1  THEN  16095 
16097  ZNI1I7I  =  U<1) 

16100  VTAB  6:  CALL  -  ’5B:  PP1H1  "FREeUENCY  SWEEP  ‘Y/Hl  ’  GOSUB  10000:ZN(56)  =  1:  IF  ZEXIOl  •  1  THEN 
16100 

ivjiA  IV  f?iji  :  j  ;»!’><!  i  jNiti:t|  -  GmSuB  JU'WiZP'iOl  '  ZfHlliil  =  15:  tUiUI*  .MOW:  t>OiU 
16150 

16110  VTAB  »:  PRINT  "FRECUENCY  M1NUHUM  : ZZ «3)  *  ZH!35l:ZZ(4l  *  IE  -  5:.’Z<5>  *  65535:ZZI6>  *  0:  GOSUB 
10100:  IF  ZEXIOl  =  1  THEN  161 10 
16120  ZNI35I  =  ZZi 31 : (1  *  35:  GOSUB  17000 

16130  VTAB  12:  PRINT  ‘FREOUENCT  HAIINW  ‘y.llO  -  ZMI34l:ZZ(4l  ;  ZHI75I-.ZZI51  *  65535: ZZ (61  *  0:  GOSUB 
10100:  IF  ZEXIOl  ?  1  THEN  16130 
16140  ZHi3«>  =  ZZ ( 3 1 : 11  -  34:  GOSUB  17000 

16150  ZNi55i  --  I:  IF  ZH(34i  >  32767  THEN  ZN'25'  =  !:  GOTO  16200:  REN  NO  HARMONIC  POSSIBLE 

ytAA  v.  fiii  .  tji.  p&tNI  ‘upMnNir  gtfcp  ■>  UKijb  j'Viiyi.  if  jcj «a>  =  j  THEN  !4!50 

16160  IF  ZZlH  -  I  THEN  ZMII22!  ■  ZH(25):!M(!23l  =  ZNI25':  6QTQ  16200 


14170  VTAB  9:  PRINT  'INITIAL  HARNOHIC  I  ';:ZZ(3l  =  ZH(l22l:ZZ(4l  =  1:ZZ(5)  =  INT  (45535  /  ZN(34)):ZZ( 
41-1:  60SU6  10100:  IF  ZEX(O)  =  1  THEN  14170 
14160  ZM( 1 22)  =  ZZI3I 

14190  VTAB  12:  PRINT  'FINAL  HARNONIC  •  '; : ZZ43)  =  ZN(123)iZZ(4)  =  ZH <  1 22) ; ZZ 45>  --  INT  (45535  /  ZNI34) 
):ZZ(41  =  1:  60SUB  10100:  IF  ZEKOI  =  1  THEN  14190 
14195  ZN4123)  =  ZZ(3):ZN(55)  =  ZHU23I  -  ZN(I22)  +  1 

14200  VTAB  4:  CALL  -  95B:  PRINT  'TINEO  SNEEP  (Y/N)  ’  ';:  GOSUB  10000:  IF  ZEX(O)  =  1  THEN  14200 
14210  IF  ZZI31  =  I  THEN  ZN(I24)  =  0:ZN<125)  =  ZH(I241:ZH<124>  =  ZN(l25»:  GOTO  14290 
14220  GOSUB  1B000:  VTAB  9:  PRINT  ZD«<0>:  REN  CLOCK 

14230  VTAB  12:  PRINT  'START  NEASURENENTS  i  HH.NN  *::ZZ(3)  =  ZN (124>:ZZ(4>  =  0:ZZ(5I  -  24:ZZ(4I  -  0:  GOSUB 
10100:  IF  ZEKOI  =  1  THEN  14230 
14240  ZNU24I  --  ZZ(3l 

14250  VTAB  15:  PRINT  'STOP  NEASURENENTS  »  HH.NN  '; : ZZ43>  -  ZHI125): ZZ(4)  =  ZN(124I:ZZ(5)  =  99.59:ZZ<4> 

=  0:  GOSUB  10100:  IF  ZEKOI  =  1  THEN  14250 
14240  ZNI125I  -  ZZ(3I 

14270  VTAB  IB:  PRINT  'HEASURENENT  INTERVAL  HH.NN  '; :ZZ<3)  -  ZN(124liZZ(4)  =  Of ZZ(5>  =  24:ZZI4>  =  0:  GOSUB 
10100:  IF  ZEKOI  =  1  THEN  14270 
142B0  ZN U 26 )  =  ZZI3I 

14290  7NI54I  =  INT  <Z«  I  (ZN(55)  I  ZH457)  1  1  :  GOTO  15300 
14300  HONE  :  VTAB  3:  PRINT  'SCALE  FACTORS:' 

14310  VTAB  4:  PRINT  'll  RESISTANCE  STANOARO  : ZZ<3>  =  ZN(51):ZZ(4)  *  .1:ZZ(5)  *  1E9:ZZ(4I  =  0:  GOSUB 
10100:  IF  ZEKOI  =  1  THEN  14310 
14320  ZN(5I)  =  ZZ(3) 

14330  VTAB  10:  PRINT  '//  CAPACITANCE  STANOARO  ';:ZZ(3I  =  ZN(52l:ZZ14l  =  0:ZZ(5)  *  IE3:ZZ(4)  -  0:  GOSUB 
10100:  IF  ZEKOI  --  I  THEN  14330 
14340  ZN452)  =  ZZ(3I 

14350  VTAB  14:  PRINT  'INPEOANCE  SCALE  FACTOR  *j:ZZ(3)  =  ZN453I :ZZ(4)  -  1:ZZ<5>  -  1E9:ZZ<4»  -  0:  GOSUB 
10100:  IF  ZEKOI  =  1  THEN  14350 
14340  ZN(53)  =  ZZ(3> 

14490  GOTO  15300 

14500  HONE  :  VTAB  3:  PRINT  'HEASURENENT  NODE:' 

14510  VTAB  5:  CALL  -  95B:  PRINT  'INPUT(S):  0=2/1':  NTAB  14:  PRINT  *1=1':  HTAB  14:  PRINT  '2=2':  HTAB 
14:  PRINT  ' 3=12-11/1':  HTAB  14:  PRINT  '4-1/ <2-11 ':  HTAB  14:  PRINT  '5=112* 

14520  VTAB  11:ZZ<3>  =  ZN <  1 27 1 : ZZ <4 1  =  0: ZZ(5>  =  5:ZZ(6I  =  1:  GOSUB  10100:  IF  ZEKOI  =  1  THEN  14520 
14525  ZN ( 1 27)  =  ZZI3I 

14530  VTAB  13:  PRINT  'OPERATOR:  0=UNITY':  HTAB  14:  PRINT  'l=JN':  HTAB  14:  PRINT  *2=I/JN':  HTAB  14:  PRINT 
'3=JN*2':  HTAB  14:  PRINT  '4=1/JN'2' 

14540  VTAB  19:?Z(3)  =  ZN(1I4):ZZ(4)  =  0: ZZ 45)  =  4:2Z<4>  =  0:  GOSUB  10100:  IF  ZEKOI  =  1  THEN  14540 
14545  ZNI114)  =  ZZ 13) 

14550  VTAB  20:  PRINT  'COORDINATES  0=A/B':  HTAB  14:  PRINT  ' 1 =R/THETA* :  HTAB  14:  PRINT  '2=L0G.R/TH£TA' 

14540  VTAB  Z3:ZZI3>  =  ZN< 46>: ZZ (4)  =  0:ZZl5)  =  2:ZZ(4l  =  0:  GOSUB  10100:  IF  ZEKOI  =  1  THEN  14540 

’H'U 

14570  GOTO  15300 

14400  HONE  :  V1AB  3:  PRINT  ’RETRIEVE  SET-UP  FILE- 
14420  GOSUB  144B0:  ZZ(B>  *  77:  GOSUB  11300:  GOTO  I  WOO 
14440  HONE  :  VTAB  3:  PRINT  'SAVE  SET-UP  FILE  TO  DISK' 

16640  GOSUB  144B0: ZZIBl  *  77:  SOSUB  11400:  GOTO  15000 

16480  1  =  FRE  10):  GOSUB  11500:  PRINT  :  PRINT  'FILE  NAHE  BILL  BE  PRECEDEO  BY  OHSSlPl-':  PRINT  :  INPUT 
‘ENTER  FILE  NANE  *  ONSS (P 1 : ZO* (01 : ZF*  =  'ONSSlPI-'  ♦  ZO*(OI:ZO*(0)  =  ZF»:  RETURN 
16-'00  HONE  :  VTAB  3:  PRINT  'EIECUTE  SET-UP  flLE' 
io; iO  rGR  ii  -  1  10  4:  GOSUB  17000:  NElT  II 


16720  II  =  M:  BOSUB  17000:  IF  ZNU4)  =  1  THEN  FOR  II  =  11  TU  13:  BOSUB  17000:  NEXT  II 
16730  II  =  ?0:  BOSUB  16900 

16740  II  -  21:  BOSUB  17000:11  *  24:  IF  ZNI23)  >  0  THEN  II  =  23 
I 6 750  BOSUB  17000:11  =  25:  60SUB  17000:11  =  26:  60SUB  17000 
16760  FOR  II  •  27  TO  30:  FOR  JZN  *  I  TO  2:  BOSUB  17100:  NETT  JZN:  NEXT  II 
16770  II  *  34:  BOSUB  17000:11  *  35:  60SUB  17000:11  =  46:  BOSUB  17000 

16780  ZNI45)  =  0201:  IF  ZNI127)  >  0  THEN  ZHM5I  =  0100:  IF  ZHII27)  >  1  THEN  ZNI451  =  0200:  IF  ZHI127)  > 

2  THEN  ZN (45)  =  0100 

16790  II  =  45:  BOSUB  17000:11  =  6B:  605UB  16900 

16B00  JZN  *  0:11  =  B5:ZH(B5I  =  1:  BOSUB  17200:  REN  DISPLAY 

16B02  JZN  *  1:ZN(B5)  =  0:  BOSUB  17200 

16B04  JZN  =  2: ZN(B5>  *  1:  BOSUB  17200 

16B06  JZN  ;  3: ZN(B5)  =  1:  60SUB  17200 

16B10  JZN  =  1:11  =  B6:ZN(B6)  -  0:  BOSUB  172O0:  REN  DISPLAY 

16B11  JZN  *  2; ZNIB6)  -  0:  BOSUB  17200 

16B12  JZN  =  3: ZN <B6)  =  Z N  < 1 16) :  BOSUB  172^0 

16B13  JZN  =  4:ZN(B6)  -  0:  BOSUB  17200 

16B14  JZN  =  5: ZN 4B6)  =  0:  BOSUB  17200 

16B16  II  :  10B:  BOSUB  17000:11  =  109:  BOSUB  17000 

16B20  60T0  1000 

16900  RESTORE  :  FOR  I  *  1  TO  II:  READ  HI;  NETT  I:Nt  =  's'*  ♦  HI  ♦  ■'*:  BOSUB  21000:  RETURN 

17000  RESTORE  :  r0R  I  *  I  TO  II:  READ  Nl:  NETT  I:NI  *  ■s’*  ♦  N»  ♦  STRI  (ZNIIII)  ♦  BOSUB  21000:  RETURN 

17100  FOR  JZN  -  1  TO  2:  RESTORE  :  FOR  I  -  1  TO  II:  READ  HI:  NETT  l:Nt  *  V*  +  Nl  ♦  STRI  (JZN)  ♦  V  ♦ 

STRI  (ZHIID)  ♦  BOSUB  21000:  NETT  JZN:  RETURN 

17200  RESTORE  :  FOR  I  =  1  TO  II:  READ  Nt:  NETT  l:Nt  =  V*  ♦  HI  ♦  STRI  (JZN)  ♦  \*  ♦  STRI  (ZN(ID)  + 

60SUB  21000:  RETURN 
IBOOO  REN  CLOCK 

1B010  ZZ(3)  -  PEEK  (  -  163B4)  -  I2B:  REN  KEY  PUSN 

19020  IF  ZZ(3)  -  1  THEN  BOTO  1000:  REN  ABORT 

180B0  IF  Z Z  (3)  =  19  THEN  IB220:  REN  STATUS 

1B090  IF  IH 3)  =  20  THEN  ZH(47)  =  999:  REN  TERNINATE 

1B100  IF  IH 3)  =  21  THEN  ZH(47)  *  0 

1B190  RETURN 

1B200  REN  TINE 

1B220  II  s  i;Hl  =  V’THO":  60SUB  24000:ZZ<3)  *  ZZ(II):NI  *  '?TN1* ’:  605UB  24000:ZZ(4)  *  ZZ(II):ZDI( 

0)  *  ‘TINE  IS  1  ♦  STRI  ( Z Z (3) >  ♦  ViNI  •  STRI  (ZZ(4)):  IF  ZZ(4)  <  9  THEN  Nl  *  T  ♦  Nl 
1B240  ZDKO)  *  ZOI(O)  ♦  HI:  VTAB  1:  HTAB  2B:  PRINT  ZDKO):  RETURN  :  REH  READ  CLOCK 
19900  REH  CHAIN  SUBROUTINE 

19920  HOHE  :  VTAB  9:NPI  *  1:  PRINT  ‘LOADING  DNSS.’ZDKO) 

19940  PRINT  CHRI  (41 i’BLOAD  CHAIN, A520,S’SBX:",D1,V’VBI 
19960  RETURN 
20000  III  *  12 

20010  ZDK4)  =  mV  ♦  CHRI  (95)  ♦  CHRI  (ZZI  ♦  64):ZDK5)  =  *J‘  ♦  CHRI  (95)  ♦  ■?*  ♦  CHRI  (ZZI  ♦  32)  ♦ 

•i  1 

20020  ZDI (6)  =  CHRI  (4)  ♦  •1NI0,:ZDI(7)  *  CHRI  (4)  ♦  'PRIO* 

20030  ZDKBI  *  CHRI  (4)  ♦  ,PRI3,:AAI  *  V  ♦  CHRI  (95)  ♦  Ms’^BBI  -  ,I'3‘>:aH: ‘ 

20900  RETURN 

21000  REH  SEND  DATA  TO  FRA 
21030  PRINT  ZDKB);  PRINT  ZDII5) 

21040  PRINT  Nl 


21470  PRINT  ZD!I7) 

21S00  RETURN 

24000  REN  DATA  RETRIEVE 

24930  PRINT  ZD! IB):  PRINT  ZD! (5) 

24040  PRINT  CHR!  113):  PRINT  H$:  6DS*'B  2BS00:  PRINT  ZD!(7);  PRINT  ZD!I6):  RETURN 

28500  PRINT  CHR!  (I3);ZD!I4);V:  PRINT  ZD4I7):  VTA8  23 

28501  ON  II  -  1  6DTD  28504,28506,28508 

28502  INPUT  ZZ(1):  PRINT  ZD!<8):  RETURN 
28504  INPUT  ZZ(I),ZZ(2):  PRINT  2D!(6):  RETURN 
28506  INPUT  ZZII),ZZI2),ZZI3):  PRINT  ZM(8)<  RETURN 
28508  INPUT  ZZI1) ,ZZ<2),ZZ(3I,ZZ(4):  PRINT  ZD4I8):  RETURN 
29000  REM  KEITHLEY  NODfcL  192  SET-UP 

29200  ZDtIO)  =  M!  +  *R0K0Q0S3N0Z0N1* 

29220  VTA8  23 

29240  PRINT  ZD!I8):  PRINT  AA!{ZH(0);8I!:  PRINT  ZD!I7):  INPUT  ZD!(0):  PRINT  ZDM);  VTA8  23:  PRINT  *1 
29260  ZZ!3)  «  VAL  I  Mil!  IZMI0),5)):  RETURN 


Appendix  B 


DIFFUSION  IN  AN  ANGULAR  CRACK 


B  -I 


This  appendix  gives  the  solution  concentration  as  a  function  of  time 
and  position  in  an  angular  crack.  The  mathematical  representation  of  the 
crack  is  shown  in  Figure  B.l.  A  crack  tip  rupture  and  dissolution 
provides  the  concentration  C  *  1  up  to  the  r  ■  a  cell  element  of  the 
angular  crack.  In  this  case,  rm  is  the  distance  from  the  crack  origin  to 


the  mc^  cell  segment  and  A  is  the  distance  between  cell  segments.  Thus, 


r  *  mA 
m 


(B-1) 


We  need  to  solve  for  cm(t)  where  cm(t)  =  average  concentration  in  rm_j  < 


r  <  r_  at  time  t  and 


in  *  1,2,3... 


Let  c  *  c(a,r,t)  be  the  solution  of  the  standard  diffusion  equation 


Sc 


n  i  a  ,  Sc. 

St  D  rSr  (r  Sr5 


(  B— 2  ) 


where  D  ■  diffusion  coefficient  with 


1  for  r  <  a 


0  for  a  <  r 


at 


Assuming  only  nearest  neighbors  interact  over  short  time  t  gtvea 


Cl(t)  -  c2(0)  +  [c L (0)  -  c2(0))  J_(l) 


(B— 3) 


and  for  m  -  2,3,4,...;  the  solution  for  the  concentration  aa  a  function  of 
time  and  position  along  the  crack  is  given  by 


B-2 


c»(t>  -  C«+1(0>  +  [c»(0)  •  W0)1  J-(n)  +  [C.-1!°>  -  c«(0»\  (n-1> 


(B-4) 


where  J_(m),  a  function  of  t  as  well  aa  m,  is  the  average  concentration  in 
cell  m  just  to  the  left  of  a  unit  jump  in  concentration  at  r  ■  rm  and  J+(m)  ia 
the  average  concentration  in  cell  (m+1)  juat  to  the  right  of  the  unit  jump  at 
r  -  rm.  That  is. 


rnA 


/  c(mA,r ,t)rdr 

J_(m)  _  - 


(B-5) 


1UU 

■ff  m— 


rdr 


(m-l)A 


-  * 

I  c(mA,r ,c)  rdr 

t  fm.  _J(°H-1)A _ 

J+'  ^  (m+1) A 


(B-6) 


JmA 


rdr 


Approximate  valuea  valid  for  /dF/A  <  <  1  are 


T  ,  .  ,  ^Dt  2m 

J_(m)  ■!  j— 


(2m  -  l)/« 


J+(n> ' 


2m 


(2m  +  l)/n 

The  derivation  of  solution  (B-4)  ia  given  aa  follows* 


Let  c  ■  Laplace  transform  of  c;  therefore. 


(  B  —7  ) 


(B-8) 


/  e  pt  c  dt 

"A 


(B-9) 


Then  the  transform  of  the  basic  problem  is 


(B-10) 


where  cQ0,  a  function  of  r.  Is  the  Initial  value  of  c: 

_  j  1  r  <  s 
coo  “  |0  a  <  r 

Solving  for  r  <  a  yields 

c  -  AIQ(qr)  +  BKQ(qr)  +  1/p  (B  -ll) 

where  IQ  and  KQ  are  the  modified  Bessel  functions;  q  ■  / p/D;  and  A  and  B  are 
constants  of  Integration*  Now  B  •*  0  because  KQ  ®  as  r  +  0. 

Solving  for  a  <  r  gives 

c  -  HIQ(qr)  +  GKQ(qr)  (B-12) 

where  H  and  G  are  constants  of  Integration  and  H  -  0  because  I  •  as 

r  + 

Both  c  snd  dc/dr  are  contlnuuus  at  r  ■  a.  Therefore, 

AlQ(qa)  +  1/p  -  GKQ(qa)  (B-13) 

Al'(qs)  -  CKo(qa)  (B-14) 

Then 

B-5 


3 

*5 


A 


qaKo  (qa) 


qa  Iq  (qa) 


by  use  of  K  (x)I  -  I  (x)K  (x) 
o  o  o  o 


1/x.  Therefore, 


l  K1(qa)  IQ(qr) 


r  <  a 


qa  1^ (qa)  KQ(qa) 


a  <  r 


Early  time  approximation  la  as  follows: 


V'Dt 


TT  <  <  l~*T<  <  1 


qa 


<  <  1 


Then 


Ko(qa)  i  Kl(qa)  -  J-fc 


Io(qa)  i  Ij (qa) 


,qfl 


/  2xqa 


1_  1_  £  -q(a-r) 

,p  2p  nJ  r  G 


r  <  a 


1  /a  -q(r-a) 

^  J7  e  a  <T 


B-6 


( B-l 5 ) 
(B-16) 


(B-17) 


(B-18) 


(B-19) 


(B-20) 


(B-21) 


(B-22 ) 


Inversion  of  the  transforms  gives 


(B-23) 


(B-24) 


With  J+(m)  and  J_(m)  defined  after  equation  (B-4)  and  with  I+(m)  and  I_(m) 
being  the  Integrals  of  the  concentration  In  those  definitions,  one  finds 


where  a  ■  mi.  The  order  of  accuracy  la  l/Q.  Approximations  correct  to  this 
order  can  be  found  as  follows: 


J  (m)  -  l - - - ■-  I  (m) 

(2m  -  1)A4 

Q&  .  r- 

I_(m)  '  f  T  IT  erfc  - - -  rdr  (B-26) 

J(  m-l)A  ^  2/  Dt 

Let  Q  ■  A/2/Dt.  Then  Q  >  >  1  by  assumption.  Integrations  are  needed  to 
order  of  accuracy  l/Q. 


Function  lerfc(x)  is  defined  by 


r  1  -x 

ierfc(x)  »  /  erfc(x)  dx  ■  - e  -  xerl'c(x) 

'tr  /iT 


a+A 


I+(m)  -  f  y  /ar  erfc  — - * 

Jr 


dr  a  “  mA 


2/Dt 


Let  x  *  (r  -  a)/2/Dt.  Then  r  ■  a  (1  +  x/mQ)  and 


I+(m)  -  f  ~  a  /I  +  x/mQ  erfc(x)  dx 


2  .Q 


2mQ 


/l  +  x/mQ  erfc(x)  dx 

-  rt 


mA 


2Q 


(H+(m)  +  K+(m)] 


where 


.1 


H  -j  / 1  +  x/mQ  erfc(x)  dx 
J0 
2 


K+  ■  (  /l  +  x/mQ  erfc(x)  dx 


H+  can  be  approximated  by 


H+  "1  +  7  mQ  +  ***  1  erfc(x> 

-  o  ^ 

*1  erfc(x)  dx 


dx 


B-8 


(B-2?) 


(B-28) 


where 


a  -  mA 
Q  -  A/2/DtT 


Then,  with  z  -  (a  -  r)/2  /Dr,  r  -  a(l  -  x/raQ).  Thus, 


2  Q 

[_(m)  -  |  /l  -  x/mQ  erfc(x)  dx 


I  (m)  - 


a  1 


2mQ 

to  lowest  order,  as  before« 


J  (m)  - 


m 


(2m  •-  l)Q/x 


(B-32) 


(B-33) 


B-10 


. *  *  *  ,•  * 


Appendix  C 


DC  TRANSMISSION  LINE  MODEL 


1  REM  "SCC  TLDC  16/ April/1 985 

2  REM  "Transmission  Line  Model  for  Stress  Corrosion  Cracking 

3  REM  "BY  M.C.H.McKubre 

10  PRINT  CHR$  (13);  CHR$  (4);"SAVE  SCC  TL  DC#60" 

500  DIM  Zn(9,200)Zm(26) 

520  F  -  96478 

600  D  ATA14..IDC,  AMPS, 0,W1DTH,CM,10,.HE1GHT,CM, 5, .LENGTH, CM,. 1, 

.ANGLE,RAD.,.02,.BULK.pH,~,6,.RHO-E,OHM-CM,20,.RKO-S,OHM-CM,lE-5, 

.G  AMM  A, CM,  1E4,#.EV  ATOM, 2, INTER  VAL.SEC.,3600, NUMBER, SEGMENTS, 20, 

TOLER  ANCE,MV,.1,MAX.1TER., LOOPS, 99 
620  DATA  26,IO.M/MftA+,A.CM-2,lE-3,V..M/MAA+,VOLT,0,lO.M/MO,A 

.CM-2,1E-5,V..M/M0,V0LT,.1,10.H/H2,A.CM-2,1E-5,V..H/H2,V0LT,.2, 

K,-— ,lE-6, MONOLAYER, COUUCM2,lE4(N*F*MW*RHO,COUUCMA3,8.4E7,D, 

CM2.S- 1 ,2E-5,PASS1  VE.l,  A.CM-2,  lE-6,?,?.0 
700  RESTORE  :  READ  II:  FOR  1  -  1  TO  11:  READ  ZQSZTS,Zm(l):  NEXT  1: 

READ  U:  FOR  1  =■  1  +  II  TO  U:  READ  ZQ$ZT$Zm(l):  NEXT  1:N%  -  Zm(12) 

1000  TEXT  :  HOME  :  VTAB  3:  PRINT  "PROGRAM  OPTIONS  ARE::  PRINT 
1010  PRINT  :  PRINT  T....SET  UP  NEW  CONSTANTS" 

1020  PRINT  :  PRINT  "2... .CALCULATE  STEADY  STATE" 

1030  PRINT  :  PRINT  "3....1N1TIALGUESS  OF  V" 

1040  PRINT  :  PRLNT  "4.  ..SAVE  ARRAY  TO  DISK" 

1050  PRINT  :  PRINT  "5....SEND  ARAY  TO  PIPELINE" 

1090  PRINT  :  PRINT  :  PRINT  "WELL  ? ";:  GET  ZQ$:  IF  ZQS  -  "0"  THEN  HOME  : 

VTAB  15:  PRINT  "PROGRAM  EXITED  -  GOTO  1000  TO  RESUME":  END 
1 100  IF  VAl  (2jQ$)  >  5  THEN  1000 
1120  ON  VAL  (ZQS)  GOTO  1500,2000,3000,4000,5000 
1140  GOTO  1000 
1500  GOSUB  1510:  GOTO  1000 

1510  HOME  :  VTAB  2:  PRINT  "SET-UP  CONSTANTS  -  PAGE#1" 

1520  RESTORE  :  READ  11:  FOR  1  -  1  TO  11:  READ  ZQSZTSJ:  VTAB  (3  +  1):  HTAB  4: 

PRINT  1;"  “ZQS;:  HTAB  20:  PRINT Zm(l);:  HTAB  PRINT  ZT$:  NEXT  1 
1540  VTAB  23:  INPUT  "HIT  #  TO  MODIFY  ";ZQS: 

IF  VAL  (ZQS)  <  1  OR  VAL  (ZQS)  >  D  THEN  1600 
1560  LI-  VAL  (ZQS):  FLASH  :  VTAB  01  +  3):  HTAB  (1):  PRINT  VTAB  22:  RESTORE  : 

READ  N:  FOR  1  -  1  TO  11:  READ  ZQS,ZTS,N:  NEXT  1:  NORMAL :  PRINT 
"CHANGE  ";ZQS;:  INPUT "  TO  "ZQS:  IF  ZQS  <  >  ""  THEN  Zm(II)  -  VAL  (ZQS) 

1580  VTAB  22:  HTAB  32:  PRINT  ZTS:  GOTO  1500 

1600  HOME  :  VTAB  2:  PRINT  "SET-UP  CONSTANTS  -  PAGE#2" 

1620  RESTORE  :  READ  II:  FOR  1  -  1  TO  11:  READ  ZQSZTSJ:  NEXT  I:  READ  U: 

FOR  1  -  11  +  1  TO  U:  READ  ZQS, ZTS, N:  VTAB  (3  +  1  - 11):  HTAB  4:  PRINT  1;" 

;ZQS;:  HTAB  20:  PRINT  Zm(l);:  HTAB  32:  PRINT  ZTS:  NEXT  1 
1640  VTAB  23:  INPUT  "HIT  #  TO  MODIFY  "ZQS:  IF  VAL  (ZQS)  -  0  THEN  RETURN 
1650  IF  VAL  (ZQS)  <  11  OR  VAL  (ZQS)  >  U  THEN  1510 

1660  U  -  VAL  (ZQS):  VTAB  22:  RESTORE  :  READ  N:  FOR  1  -  1  TO  11:  READ  ZQS.ZTS.N:  NEXT  1: 
READ  N:  FOR  1  -  II  +  1  TO  U  READ  ZQSZTS.N:  NEXT  1:  NORMAL  :  PRINT  "CHANGE 
"ZQS;:  INPUT  "  TO  "ZQS:  IF  ZQS  <  >  THEN  Zm(U)  -  VAL  (ZQS) 

1680  VTAB  22:  HTAB  32:  PRINT  ZTS:  GOTO  1600 

2000  HOME  :  VTAB  2:  PRINT  "CALCULATING  STEADY  STATE  CONDITIONS’ 

2010  PRINT  CHRS  (13);  CHRS  (4);"PR#r:  PRINT  CHRS  (15): 

PRINT  CHRS  (9);"132N":  PRINT  CHRS  (12) 

2012  RESTORE  :  READ  11:  FOR  1  -  I  TO  11:  READ  ZQSZTSJ:  PRINT  1;")  "ZQS;"  -  "Zm(l):"  "ZTS: 
NEXT  1:  READ  U:  FOR  1  -  II  +  1  TO  U:  READ ZQSZTS.J:  PRINT  1;")  "ZQS;*  -  ";Zm(l);"  "ZTS: 
NEXT  1:  PRINT  :  PRINT  "VERSION: ";:  LIST  10 

2014  ZQS  -  "  /  ":  PRINT  "LOOP" ZQS;*NQ"ZQS;"G"ZQS;"SI"ZQS;*rZQS;“V"ZQ$  -  "  " 


02 


APPENDIX  ??  DC  Transmission  Line  Model 


2020  :SI  =  Zm(l) 

2040  N%  -  Zm(12):  REM  "#  OF  SEGMENTS 
2044  C  -77.4:  REM  "F/.5RT 
2046  T  -  0:  REM  "TIME-0 

2100  Zn(0,0)  -  0:Zn(7,0)  -  10  A(-  Zm(6)):Zn(8,0)  -  Zn(7,0)  AZm(I0)  /  Zm(21):  FOR  N  -  I  TO  N%: 

FOR  M  =  0  TO  6:Zn(M,N)  -  0:  NEXTM:Zn(7,N)  -  Zn(7,0):Zn(8,N)  -  Zn(8,0):Zn(9N)  =  I: 

NEXT  N:  REM  "INITIALISE  ARRAY 

2120  IF  Zn(l,0)  =  0  THEN  Zn(l,0)  =  (Zm(I6)  +  Zm(20))  /  2:  REM  "INITIAL  GUESS 
2200  LL  -  0:G  -  32:NQ  -  0:  REM  "LOOP  COUNTER  &  INITIAL  VSTEP  (MV) 

2220  SI  -  Zm(l):I  =  0:LL  -  LL  +  I:  FOR  N  -  I  TO  N% 

2240  GOSUB  10000:TI  -  Zn(4,N)  +  Zn(5,N)  +  Zn(6,N): 

REM  "CALCULATE  INTERFACE ALCURRENTS 
2250  Zn(2,N)  -  Zm(8)  *  (Zm(4)  /  N%)  /  (Zm(2)  *  Zm(3)):  REM  "RS 

2260  Zn(3,N)  -  Zm(7)  /  (2  ♦  Zm(2)  ♦  (N%  +  1  -  N)  *  TAN  (Zm(5))):  REM  "RE 

2270  Zn(0,N)  =  Zn(0,N  -  1)  +  Zn(2,N)  ♦  (SI  -  IT):  REM  "VS 

2280  Zn(l,N)  -  Zn(l,N  -  1)  -  Zn(3,N)  ♦  (SI  -  TI):  REM  ”VE 

2290  I  - 1  +  TI:  REM  "MOVING  SUM  I 

2299  NEXT  N 

2300  IF  ABS  (G)  <  Zm(  13)  THEN  2500:  REM  "SUCCESSFUL  ITERATION 

2301  ZQ$  = "  ":  PRINT  UL^Q$;NQ^QS;G2Q$;SI2Q$;I2Q$(Zn(l,0) 

2305  GOSUB  9000:  REM  "READ  KEYBOARD 

2310  IF  I  <  SI  THEN  2400 

2320  REM  "VOLTAGE  TOO  LARGE 

2330  IF  NQ  -  1  THEN  G  -  G  /  2 

2340  NQ  -  - 1:  REM  "LAST  LOOP  FLAG 

2350  Zn(I,0)  -  Zn(I,0)  -  G  /  1000:  GOTO  2220 

2400  REM  "VOLTAGE  TOO  SMALL 

2410  IF  NQ  *  - 1  THEN  G  -  G  #  2 

2420  NQ  -  l:Zn(l,0)  -  Zn(l,0)  +  G  /  1000:  GOTO  2220 

2500  REM  "SUCCESSFUL  ITERATION  TO  FIND  DC  CONDITION  @  TIME  T 

2520  PRINT  :  PRINT  "  V  -  ";Zn(l,0);"  VOLTS,  I/I.SUM  -  ";Zm(l AMPS";",  ERROR  - 

";G;"MV,  TIME  -  ";LT;"  HOURS,  LENGTH  -  ";Zm(4);"  CM":  PRINT  "  CRACK  GROWTH  - 
";DL/Zm(ll);"  CM/S,  TIP  PH-";  -  LOG  (Zn(7,N%))  /  LOG  (10):  PRINT 
2540  LT  -  LT  +  Zm(ll)  /  3600:  REM  "TIME 

2549  PRINT  "BEFORE  DELTA-C":ZQS  - "  ":  FOR  N  -  0  TO  N%  STEP  N%  /  5: 

FOR  M  -  0  TO  6:PRINT  Zn(M,N);ZQS;:  NEXT  M:  PRINT  N:  NEXT  N:  PRINT  :  PRINT 
2560  FOR  N  -  1  TO  Nft:  GOSUB  1 1000:  NEXT  N:  REM  "NOW  SOLVE 
FOR  DELTA-CONCENTRATIONS 

2580  Zm(4)  -  Zm(4)  +  DL:  FOR  N  -  1  TO  N%:  FOR  M  -  7  TO  8:Zn(M,N)  -  Zn(M,N)  *  ((Zm(4)  -  DL) 

/  Zm(4))  A2:  NEXT  M:Zn(9,N)  -  Zn(9,N)  *  (Zm(4)  -  DL)  /  Zm(4):  NEXT  N:  REM  "C 
SPREAD-OUT  WITH  CRACK  GROW1H 

2599  PRINT  "BEFORE  DIFFUSION" :ZQS  -  "  ":  FOR  N  -0  TO  N%  STEP  N%  / 
5:PRlNTZn(0,N);ZQS-,Zn(l,N);ZQ$;:  FOR  M  -4  TO  9:  PRINT  Zn(M,N);ZQ$;:  NEXT  M:  PRINT 
N:  NEXT  N: PRINT  :  PRINT 

2600  REM  "CALCULATE  DIFFUSION 

2620  ZZ(0)  -  (N%  /  Zm(4))  *  (Zm(24)  *  Zm(ll)  *  3.14)  \5:  REM  CHI 

2640  FOR  N  -  N%  TO  1  STEP  -  1:  GOSUB  12000  NEXT  N 

2699  PRINT  "AFTER  DIFFUS10N":ZQ$  -  "  ":  FOR  N  -  0  TO  N%:  FOR  M  -  I  TO  3  STEP  2:  PRINT 
Zn(M,N)2QS;:  NEXT  M:  FOR  M  -  4  TO  9:  PRINT  Zn(M,N);ZQS;:  NEXT  M:  PRINT  N:  NEXT 
N:  PRINT :  PRINT 

2920  LL  -  0:G  -  4:NQ  -  0:  GOTO  22 

3000  HOME  :  VTAB  4:  PRINT  "INITIAL  GUESS  OF  VDC  -  ";Zn(l,0):  PRINT  :  INPUT  "CHANGE 
T0"2QS:  IF  ZQS  -  "0"  THEN  1000 

3020  IF  VAL  (ZQS)  <  >  0  AND  ABS  (  VaL  (ZQS))  <  2  THEN  Zn(l,0)  -  VAL  (ZQS) 

3040  GOTO  1000 
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5000  REM  "OUTPUT  DATA  TO  PIPELINE 
5010  GOSUB  5020:  GOTO  1000 

5020  PRINT  CHR$  (13):  CHR$  (4);"PR#1";  PRINT  CHR$  (9);"132N":  PRINT  CHR$  (15) 

5040  FOR  N  -  0  TO  N%:  FOR  M  -  0  TO  6:  PRINT  Zn(M,N);"  " NEXT  M:  PRINT  N:  NEXT  N 

5080  PRINT  CHR$  (13);  CHR$  (4);"PR#0":  RETURN 

9000  REM  "READ  KEYBOARD 

9010  ZZ(3)  »  PEEK  ( -  16384)  -  128:  POKE  -  16386,0:  IF  ZZ(3)  <  1  THEN  RETURN 

9020  PRINT  CHR$  (13);  CHR$(4);"PR#0" 

9030  IF  72X1)  -  16  THEN  PRINT "':  GOSUB  1510:  REM  "CHANGE  PARAMETERS 

9040  IF  ZZ{3)  =  20  THEN  POP  :  GOTO  1000 

9090  PRINT  CHR$  (13);  CHR$  (4);"PR#1":  IF  ZZ(3)  -  16  THEN  GOSUB  9100 

9099  RETURN 

9100  REM  "PRINT  PARAMETERS 

91 10  RESTORE  :  READ  D:  FOR  1  -  1  TO  11:  READ  ZQS^TSJ:  PRINT  1;")  "(ZQ$;"  -  ”;Ziti(1);"  ";ZTS: 
NEXT  1:  READ  U:  FOR  1  -  D  +  1  TO  U:  READ  ZQ$,ZTSJ:  PRINT  1;")  "ZQS;"  -  ";Zm(I);"  "2T$: 
NEXT  1:  PRINT  :  RETURN 

9999  PRINT  CHR$  (13);  CHR$  (4);"SAVE  SCC  TL  DC#6A":  END 

10000  REM  "INTERFACIAL  ELECTRICS 

10100  REM  "METAL  DISSOLUTION  @  CRACK  TIP  DUE  TO  STRESS 
10110  Zn(4,N) -0:  IF  N  <N%  THEN  10200 

10120  ZZ(0)  =  Zm(l5)  *  Zm(9)  *  Zm(2)  ♦  TAN  (Zm(5)):ZZ(l)  -  Zm(10):ZZ(2)  -  Zm(16):ZZ(3)  - 
Zn(8,N):  GOSUB  10800:  REM  "B-V  0-10,I-N>V0,3-C 
10140  Zn(4,N)  -  ZZ(5):DL  -  Zn(4,N)  •  Zm(l  1)  /  (Zm(23)  *  Zm(9)  *  Zm(5)  *  Zm(2)):  IF  DL  <  0  THEN 

DL  -  0:  REM  "DELTA-L 
10200  REM  "CHEMICAL  PASSIVATION 
10300  REM  "HYDROGEN  ION  REDUCTION 

10320  ZZ(0)  -  Zm(l9)  *  2  *  Zm(2)  *  (Zm(4)  /  N%):ZZ<1)  -  122(2)  -  Zm(20):ZZ(3)  -  Zn(7,N):  GOSUB 
10800:  REM  "B-V  0-I0,l-N,2-V0,3-C 
10340  Zn(5,N)  -  ZZ{5) 

10400  REM  "METAL  OXIDE  REDUCTION 
10410  ZZ(0)  -  Zn(9,N):  IF  ZZ(0)  >  1  THEN  ZZ(0)  -  1 

10420  ZZ(0)  -  ZZ(0)  *  Zm(l7)  ♦  2  ♦  Zm(2)  ♦  (Zm(4)  /  N%):ZZ(1)  -  Zm(10):ZZ<2)  -  Zm(18):ZZ(3)  - 
Zn(7,N):  GOSUB  10800:  REM  "B-V  0-10,l-N,2-V0,3-C 
10440  Zn(6,N)  -  ZZ(5) 

1 0500  REM  " METAL  DISSOLUTION  ON  CRACK  WALLS 

10510  ZZ(9)  -  ( 1  -  Zn(9,N)):  IF  ZZ(9)  <  0  THEN  ZZ(9)  -  0:  REM  "ACTIVE/PASSIVE 

10520  ZZ(0)  -  Zm(15)  •  2  •  Zm(2)  ♦  (Zm(4)  /  N%)2Z(1)  -  Zm(10):ZZa)  -  Zm(16)JC(3)  -  Zn(8, 

N):  GOSUB  10800:  REM  "B-V  0-10,1 -N.2-V0.3-C 
10540  ZZ(0)  -  ( I  -  Zn(9,N)):  IF  ZZ(0)  <  0  THEN  ZZ(0)  -  0:  REM  THETA 
10560  IF  72X0)  >  C  THEN  ZZ<6)  -  ZZ(5)  •  ZZ(0):  GOTO  10590:  REM  "ACTIVE  KINETICS 
10580  72X6)  -  SGN  (ZZ(5))  *  ABS  (Zxn(25)  *  2  ♦  Zm(2)  *  (Zm(4)  /  N%)  *  (1  -  Zn(8,N)>  /  Zn(9,N)):  IF 
ABS  (ZZ(5))  <  ABS  (2Z(6))  THEN  ZZ(6)  -  2Z(5):  REM  "PASSIVE  DIFFUSION 
10590  Zn(4,N)  -  Zn(4,N)  +  ZZ(6) 

10600  RETURN 

10800  REM  "COMPUTE  BUTLER- VOLMER  1 

10805  72X2)  -  Zn(l.N  -  1)  -  Zn(0,N  -  1)  -  72(2 )  +  (2  /  (Zm(10)  *  O)  *  LOG  (2Z<3)):  REM  "ZETA 

INCLUDING  CONC. 

10810  72X4)  -  22X1)  •  C  •  (ZZ(2)):  IF  ABS  (ZZ(4))  >  30  THEN  ZZ(4)  -  SGN  (ZZ(4))  •  30 
10820  72X5)  -  ZZ(0)  ♦  (  EXP  (ZZ(4))  -  EXP  ( -  ZZ(4))) 

10830  72X6)  -  Zm(2)  *  Zm(5)  •  (Zm(4)  /  N%)  A2  *  (2  *  N%  -  2  *  N  +  1):  REM  "VOLUME 

10839  IF  N  -  NS,  AND  ZZ(5)  >  0  THEN  10849:  REM  "NO  LIMIT  TO  M->MO  @  TIP 

10840  IF  ABS  (72X5))  >  ABS  (99  •  ZZ(1)  *  F  *  72X5)  *  ZZ(6)/  Zm(ll))  THEN  ZZ(5)  -  SGN  (ZZ<5))  * 
99  •  72XD  *  96478  •  72X5)  *  ZZ(6)  /  Zm(ll):  REM  "MAX  DELTA-C 
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10849  REM  "ZQS  - "  ":  FOR  NN  -  0  TO  6:  PRINT  ZZ(NN);ZQ$;:  NEXT  NN:  PRINT  N 

10850  RETURN 

1 1000  REM  "INTERFACIAL  CHEMISTRY 

1 1020  ZZ(0)  -  Zm(2)  ♦  Zm(5)  *  (Zm(4)  /  N%)  A2  *  (2  *  N%  -  2  *  N  +  1):  REM  "V 
11040  ZZ(1)  -  Zn(4,N)  *  Zm(ll) /  (Zm(l0)  ♦  F  *  ZZ(0)):  IF  ZZ(1)  >  Zn(8,N)  /  2  THEN  ZZ(1)  - 
Zn(8,N)  /  2:  REAf  "DELTA  METAL  IONS 

11060  720)  =  (Zn(5,N)  +  Zn(6,N))  *  Zm(ll)  /  (F  *  ZZ(0)):  IF  ZZ(2)  >  Zn(7,N)  /  2  THEN  ZZ(2)  - 
Zn(7,N)  /  2 JfEAf  "DELTA  PROTONS 

1 1080  720)  -  Zn(6,N)  *  Zm(ll)  /  (Zm(lO)  *  Zm(22)  *  Zm(2)  *  (Zm(4)  I  N%)):  IF  ZZ(3)  >  Zn(9,N)  /  2 
THEN  720)  -  Zn(8,N)  /  2:  REA/  "DELTA  METAL  OXIDE 
11100  Zn(7,N)  =  Zn(7,N)  -  ZZ(2):Zn(8N)  -  Zn(8,N)  -  ZZ(l):Zn(9,N)  -  Zn(9,N)  -  720):  REM : 
"ELECTROCHEMISTRY 

1 1 120  720)  -  (Zn(8,N)  -  720))  *  F  *  Zm(l0)  /  (Zm(23)  *  Zm(2)  *  (Zm(4)  /  N%)  *  2):  IF  ZZ<9) 
>Zn(9,N)  /  2  THEN  ZZ(9)  -  Zr(9,N)  /  2:  REAf  "DELTA  OXIDE 
1 1 140  ZZ(5)  -  Zn(9,N):  IF  ZZ(5)  >  1  THEN  ZZ(5)  -  1 :  REM  "MO  ACTIVITY 
11160  Zn(7,N)  -  (Zn(8,N)  *  Zm(21)  /  ZZ(5))  A(i  /  Zma0))Zn(8,N)  -  ZZ(5)  ‘  (Zn(7,N)  AZm(10))  / 
Zm(21)Zn(9,N)  -Zn(9,N)  -  ZZ(9):  REAf  "CHEMISTRY 
11600  RETURN 

12000  REAf 

12020  720)  -  2  •  <N%  -  N):ZZ(1)  -  (ZZ(3)  +  3)  /  (ZZ<3)  +  1):  REM  "V  N-l/N 
12040  720)  -  (720)  -  1)  /  (ZZ(3)  +  1):  IF  720)  <  0  THEN  720)  -  0:  REM  "V  N+l/N 
12080  IF  720)  >  -5  THEN  12200 
12100  REM  "SHORT  TIME 

12120  FOR  M  -  7  TO  8:Zn(M,N)  -  Zn(M,N)  *  (1  -  ZZ(0))  +  ZZ(1)  *  7Z(0)  *  (3  *  Zn(M,N)  +  Zn(M,N  - 
1))  /  4  +  ZZa)  *  72(0)  *  (3  *  Zn(M,N)  +  Zn(M,N  +  l»  /  8:  NEXT  M 
12160  RETURN 
12200  REM  "LONG  TIME 

12220  FOR  M  -  7  TO  8:Zn<M,N)  -  ZZ(1)  *  (3  *  Zn(M,N)  +  Zn(M,N  - 1))  /  8  +  ZZ(2)  *  (3  *  Zn(M,N) 
+Zn(M,N  +  1))  /  8:  NEXT  M 
12260  RETURN 
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